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ABSTRACT .  An  investigation  of  the  effect  of  shock 
wave-boundary  layer  interaction  on  aerodynamic 
heat-transfer  coefficients  is  reported.  A  liter¬ 
ature  survey  and  empirical  expressions  for  predicting 
peak  heat-transfer  coefficients  in  interaction  regions, 
developed  with  the  use  of  existing  experimental  data, 
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h 
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X 
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8  Shock  angle 
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6  Boundary-layer  thickness 
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p  Density 

Subscripts 

0  Beginning  of  pressure  rise 
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2  Shock  impingement,  immediately  behind  impinging  shock; 
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Section  1 
INTRODUCTION 


BACKGROUND  AND  OBJECTIVES 

The  problem  of  aerodynamic  heating,  which  has  developed  with  the 
advent  of  supersonic  flight,  can  be  compounded  considerably  by  various 
factors  affecting  the  flow  field.  One  of  the  most  prominent  of  these 
factors  is  shock  wave-boundary  layer  interaction.  In  addition  to 
causing  order-of-magnitude  increases  in  the  local  aerodynamic  heating 
rate,  shock  wave-boundary  layer  interaction  alters  the  pressure  dis¬ 
tribution  over  a  body  and  may  ca  e  boundary-layer  separation.  Shock 
wave-boundary  layer  interaction  will  occur  on  any  vehicle  with  control 
surfaces,  or  other  surface  protuberances  during  supersonic  and  hyper¬ 
sonic  flight. 

Extensive  work,  both  analytical  and  experimental,  has  been  performed 
in  the  study  of  shock  wave-boundary  layer  interaction.  Most  of  the  early 
work  was  concentrated  on  the  pressure  distribution  in  interaction  re¬ 
gions;  however,  more  recent  studies  have  investigated  the  heat-transfer 
aspects  of  the  problem.  As  the  first  phase  of  this  study,  a  literature 
survey  was  made  to  determine  (1)  the  current  state-of-the-art  of  analyt¬ 
ical  solutions  to  the  problem,  and  (2)  the  extent  and  findings  of 
experimental  investigations.  Results  ,f  this  survey  are  presented  in 
Section  3.  An  extensive  reference  list,  an  author  index,  and  a  bib¬ 
liography,  giving  related  reports  that  were  not  used  in  this  study,  are 
included  in  this  report. 

Because  the  nature  of  the  flow  fi  'Id  in  shock  wave-boundary  layer 
interaction  regions  is  extremely  complex,  the  analytical  treatment  of 
the  problem  presents  a  great  deal  of  difficulty.  A  few  successful 
methods  of  solution  for  laminar  boundary-layer  interaction  have  been 
developed  but  are  quite  involved.  No  analytical  solutions  for  turbulent 
interactions  are  known  to  exist.  The  major  difficulty  in  obtaining 
satisfactory  analytical  solutions  has  be  n  in  accurately  describing  the 
velocity  and  enthalpy  profiles  in  the  interaction  region,  particularly 
in  cases  in  which  separation  occurs.  Son; j  of  the  more  successful  ap¬ 
proaches  to  the  problem  are  reviewed  in  Faction  3. 

Because  there  are  complexities  involved  in  analytical  methods,  a 
need  exists  for  a  rapid  engineering  method  for  estimating  heat  transfer 
in  shock-interaction  regions.  Several  empirical  relations  have  been 
previously  developed  from  the  results  of  various  experimental  investi¬ 
gations;  however,  most  of  these  were  for  specific  configurations  and 
were  based  on  a  limited  amount  of  data.  The  objective  of  this  study  is 
to  correlate  as  much  of  the  existing  experimental  data  as  possible  to 
determine  the  feasibility  of  developing  general  engineering  methods  for 
calculating  heat  transfer  in  shock  wave-boundary  layer  interaction  regions. 
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TYPES  OK  SHOCK  WAVE-BOUNDARY  LAYER  INTERACTIONS 


There  are  several  manners  in  which  shock  wave-boundary  layer  inter¬ 
action  may  arise.  The  iollowing  three  types  of  interactions  are  treated 
in  this  study: 

1.  Planar  shock  wave,  generated  externally,  impinging  on  and 
reflecting  from  a  flat  plate  or  cylinder. 

2.  Planar  or  conical  shock  wave  generated  within  a  boundary  layer 
by  a  compression  corner. 

3.  Planar  shock  wave  generated  by  the  leading  edge  of  a  fin 
mounted  on  a  flat  plate  interacting  with  the  boundary  layer  of  the 
plate. 

Two  examples  of  the  impinging-shock  interaction  are  shown  in  Fig.  1. 
Figure  1(a)  indicates  a  shock  wave  generated  by  the  leading  edge  of  an 
aircraft  wing  impinging  upon  a  captive-flight  missile.  Other  sources 
of  impinging  shock  waves  could  be  the  aircraft  air  inlet,  the  missile 
launcher  assembly,  or  additional  underwing  ordnance. 


FIG.  1.  Impinging-Shock  Interaction  on  (a)  a  Captive-Flight 
Missile,  and  (b)  an  Inlet  Diffuser. 
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As  part  of  a  study  to  investigate  flow-interference  effects  in 
captive  flight,  this  Center  (NWC)  has  constructed  and  is  flight  testing 
an  instrumented  store  of  typical  airborne  configuration  (Ref.  1).  The 
round  is  currently  instrumented  to  measure  surface  pressure,  and  it  has 
been  used  to  obtain  pressure  distributions  at  speeds  up  to  Mach  1.9. 
Pressure  distributions  from  supersonic  flights  have  indicated  the  pres¬ 
ence  of  impinging  shock  waves  generated  by  the  launcher  assembly.  With 
the  use  of  the  results,  developed  later  in  this  report,  the  measured 
pressure  rises  in  the  shock-impingement  regions  correspond  in  some  cases 
to  increases  in  heat-transfer  coefficient  of  over  100%.  Increased 
heating  rates  and  thermal  gradients  in  the  interaction  regions  could 
create  serious  aerodynamic  heating  problems,  particularly  if  shock 
impingement  occurred  in  the  vicinity  of  the  missile  motor  or  warhead. 
Altered  pressure  distributions  in  the  shock-interaction  regions  could 
also  have  an  effect  on  separation  of  the  store  from  the  aircraft. 

A  problem  which  will  exist  in  hypersonic  air-breathing  propulsion 
systems  such  as  supersonic  combustion  ramjets  (Scramjets)  is  shown  in 
Fig.  1(b).  This  sketch  illustrates  the  impingement  upon  the  inlet  dif¬ 
fuser  of  a  shock  wave  generated  by  the  lip  of  the  air  inlet. 

Configurations  which  would  cause  compression  corner  interactions 
are  a  deflected  flap,  such  as  shown  in  Fig.  2,  or  a  body  with  a  conical 
flare. 


FIG.  2.  Compression-Comer  Interaction  Generated 
by  a  Deflected  Flap. 


The  type  of  interaction  shown  in  Fig.  3  is  present  on  any  body  with 
wings  or  fins  while  in  supersonic  flight.  The  leading  edge  of  the  fin 
generates  a  planar  shock  wave  which  interacts  with  the  boundary  lay.er 
on  the  surface  of  the  body. 
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Other  types  of  shock  wave-boundary  layer  interactions,  not  treated 
in  this  study,  are  shown  in  Fig.  4,  5,  and  6.  Figure  4  shows  an  inter¬ 
action  caused  by  a  surface  protuberance,  in  this  case  a  forward-facing 
step.  The  upstream  portion  of  this  type  of  interaction  is  similar  to 
that  of  the  previous  types  mentioned.  This  is  an  example  of  the  free- 
interaction  concept,  which  states  that  the  upstream  portion  of  the 
interaction  region  is  independent  of  the  source  of  the  disturbance. 
Several  analytical  and  experimental  investigations  of  this  type  of 
interaction  have  been  made  (Ref.  2,  3,  4,  5,  6,  7,  and  8).  The  shock- 
on-shock  interaction  (Ref.  9,  10,  and  11)  is  illustrated  in  Fig.  5. 

This  is  a  special  case  of  the  shock-impingement  interaction  in  which  an 
externally  generated  shock  wave  impinges  on  a  surface  that  has  a  leading- 
edge  shock.  The  externally  generated  shock  interacts  with  the  leading- 
edge  shock  before  impinging  on  the  boundary  layer.  An  interaction  of 
this  type  would  result  from  a  nose-generated  bow  shock  impinging  on  the 
leading  edge  of  a  wing  or  fin.  Hypersonic  boundary-layer  interaction 
(Ref.  12  and  13),  shown  in  Fig.  6,  arises  when  the  shock  wave  generated 
at  the  leading  edge  of  a  surface  is  at  a  small  enough  angle  to  interact 
w  th  the  boundary  layer. 


SHOCK  WAVE 


FIG.  4.  Shock  Wave-Boundary  Layer  Interaction 
Generated  by  a  Forward-Facing  Step. 
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FIG.  5.  Shock-On-Shock  Interaction. 


LEADING-EDGE 


FIG.  6.  Hypersonic  Shock  Wave-Boundary 
Layer  Interaction. 
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Section  2 

INTERACTION-REGION  FLOW  PATTERN 


SHOCK  IMPINGEMENT 

The  basic  features  of  the  shock-interaction  region  for  the  case  of 
a  planar  shock  wave  impinging  and  reflecting  from  a  laminar  boundary 
layer  are  shown  in  Fig.  7.  The  pressure  rise  induced  by  the  incident 
shock  wave  is  propagated  upstream  through  the  subsonic  portion  of  the 
boundary  layer.  As  a  result,  the  boundary  layer  thickens,  and  its 
momentum  decreases.  If  the  adverse  pressure  gradient  is  great  enough, 
the  skin  friction  will  be  reduced  to  2ero  and  the  boundary  layer  will 
separate.  The  region  of  separated  flow  is  defined  by  the  dividing 
streamline.  Flow  reversal  takes  place  below  the  zero  velocity  stream¬ 
line,  as  indicated  by  the  velocity  profile  (Fig.  7).  Streamline 
curvature  away  from  the  surface  in  the  region  of  the  thickening  boundary 
layer  creates  a  series  of  compr  3 ion  waves  that  merge  to  form  a  reflected 
shock  wave,  sometimes  referred  to  as  the  separation  shock.  The  incident 
shock  is  reflected  from  the  constant-pressure  separated  region  as  a 
Pra~  ’tl-Meyer  expansion  fan.  Upon  passing  through  the  incident  shock 
and  expansion  fan,  the  flow  is  turned  back  toward  the  surface.  As  the 
flow  reattaches  and  turns  parallel  to  the  plate,  a  second  reflected 
shock,  termed  the  reattachment  shock,  is  formed.  Immediately  downstream 
of  reattachment,  the  boundary-layer  thickness  reaches  a  minimum.  It  is 
in  this  region  that  the  maximum  heating  rates  occur. 

In  the  case  of  shock  impingement  on  a  turbulent  boundary  layer, 
the  length  of  the  interaction  region  is  considerably  shorter  than  that 
of  the  laminar  boundary  layer.  This  results  because  a  turbulent 
boundary  layer  has  greater  momentum  than  has  a  laminar  boundary  layer, 
and  it  can  more  successfully  sustain  an  adverse  pressure  gradient.  As 
a  result,  a  much  greater  pressure  rise  is  required  to  cause  separation 
of  a  turbulent  boundary  layer  than  separation  of  a  laminar  boundary 
layer. 


COMPRESSION  CORNER 

The  flow  field  corresponding  to  a  compression-corner  interaction 
is  shewn  in  Fig.  8.  As  in  the  case  of  the  impinging-shock  interaction, 
the  pressure  disturbance,  this  time  generated  by  the  compression  corner, 
propagates  upstream  through  the  subsonic  portion  of  the  boundary  layer. 
In  the  process,  the  boundary  layer  thickens;  a  compression  shock  is 
formed;  and  if  the  pressure  gradient  is  large  enough,  boundary-lay e->- 
separation  occurs.  If  the  turning  angle  is  large  enough,  the  compres¬ 
sion  will  take  place  in  two  phases  (Fig.  8).  In  the  first  phase. 
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FIG.  7.  Schematic  Representation  of  an  Oblique 
Shock  Wave  Impinging  on  a  Laminar  Boundary  Layer. 


FIG.  8.  Schematic  Representation  of  a  Compression- 
Corner-Generated  Interaction. 


7 


►  r  **v*,«4.  ',j 


NWC  TP  4485 


a  separation  shock  is  formed  as  a  result  of  streamline  curvature  caused 
by  separation;  in  the  second  phase,  a  reattachment  shock  is  formed  as 
the  flow  tu'rns  parallel  to  the  compression  surface.  Experiments  show 
that  pressure  and  heat-transfer  distributions  are  quite  similar  for  the 
impinging  shock  and  compression-corner  interactions. 


FIN-GENERA f ED  SHOCK  WAVE 

The  5  low  field  associated  with  the  interaction  between  the  shock 
wave  generated  by  the  leading  edge  of  an  unswept  fin  and  a  boundary 
layer  is  considerably  more  complicated  than  the  two  cases  discussed 
previously.  This  arises  from  the  fact  that  the  sweep  of  the  fin  shock 
creates  .1  three-dimensional  rather  than  a  two-dimensional  interaction. 
If,  however,  the  flow  field  is  viewed  in  a  plane  normal  to  the  direction 
of  the  shock  wave  (Fig.  9),  the  flow  field  appears  similar  to  that  of  a 
two-dimensional  interaction.  Although  the  lambda-shock  configuration 
differs  from  the  previous  interactions,  the  separation  shock,  expansion 
fan,  ard  reattachment  shock  are  stilL  present.  Pressure  and  heat- 
transfer  distributions  take  the  same  form  as  those  in  the  impinging- 
shock  and  compression-corner  interactions.  Leading-edge  sweep  alters 
these  distributions,  but  the  general  trend  remains  the  same. 


FIG.  9.  Schematic  Representation  of  a  Fin-Shock- 
Generated  Interaction. 
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PARAMETERS  REQUIRED  TO  DESCRIBE  INTERACTION 

To  adequately  describe  the  interaction,  the  following  conditions 
and  parameters  must  be  determined. 

1.  Type  of  boundary  layer 

a.  Laminar  (transition  downstream  of  interaction); 
turbulent  (transition  upstream  of  interaction) ;  or 
transitional  (transition  within  interaction  region) 

b.  Attached  or  separated 

2.  Interaction  geometry 

a.  Location  and  length  of  interaction  region 

b.  Separation  and  reattachment  points 

3.  Pressure  distribution  (primarily  peak  pressure  and  in 
separated  regions,  plateau  pressure) 

4.  Heat-transfer  distribution  (primarily  peak  heat  transfer  and 
heat  transfer  in  separated  region) 
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Section  3 
LITERATURE  SEARCH 


ANALYTICAL  STUDIES 

Because  the  nature  of  shock  wave-noundary  layer  interactions  is 
complex,  attempts  at  analytical  solutions  to  the  problem  have  achieved 
only  limited  success.  The  successful  approaches  are  limited  to  the 
case  of  two-dimensional  or  axisymmetric  laminar  b<-’>ndary  layers;  although 
some  semiempirical  work  has  been  performed  on  turbulent  interactions 
(Ref.  14).  Almost  all  of  the  approaches  utilize  momentum  integral 
methods.  The  flow  is  described  by  solving  the  boundary- layer  equations 
governing  conservation  of  mass,  momentum,  and,  if  the  flow  is  nonadiabatic, 
energy.  The  primary  difficulty  in  solving  the  boundary-layer  equations 
is  in  specifying  velocity  and  enthalpy  profiles  that  accurately  describe 
the  flow  reversal  which  takes  place  in  separated  regions. 

Several  of  the  earlier  studies  (Ref.  15,  16,  and  17)  utilized  the 
Karman-Pohlhausen  momentum-integral  method,  in  which  the  velocity  profile 
is  related  to  the  local  static-pressure  gradient.  Some  success  was 
achieved  by  Martellucci  and  Libby  (Ref.  17),  who  used  this  method  with 
quartic  velocity  profiles.  These  approaches,  however,  did  not  adequately 
describe  conditions  within  separated  regions. 

Improved  solutions  were  obtained  with  a  semiempirical  method  devel¬ 
oped  by  Crocco  and  Lees  (Ref.  18)  and  rodified  by  Click  (Ref.  19),  in 
which  the  velocity  profile  was  not  directly  related  to  the  local  static- 
pressure  gradient.  The  semiempirical  features  of  the  Crocco-Lees  method 
were  eliminated  by  Taui  (Ref.  20)  in  a  method  which,  in  addition  to  the 
aforementioned  boundary- layer  equations,  utilized  the  moment  of  momentum 
equation.  This  approach,  termed  the  two-moment  method,  gave  good  agree¬ 
ment  with  experimental  data  except  in  separated  regions.  Difficulties 
were  encountered  in  these  regions  because  of  the  use  of  quartic  velocity 
profiles.  Lees  and  Reeves  (Ref.  21)  demonstrated  that  polynomial  veloc¬ 
ity  profiles  described  by  a  single  parameter  do  not  exhibit  the  correct 
behavior  in  separated  flows.  The  velocity  profiles  in  a  separated  region 
can  be  more  accurately  described  by  polynomial  expressions  with  two  or 
more  parameters,  or  by  representing  the  velocity  profiles  above  and 
below  the  zero  velocity  streamline  (Fig.  7)  by  a  different  family  of 
profiles.  Lees  and  Reeves  used  the  later  approach,  obtaining  the  veloc¬ 
ity  and  enthalpy  profiles  for  attached  and  separated  flows  by  using  the 
upper  and  lower  branch  solutions,  respectively,  of  the  compressible  flow 
analogs  to  the  Falkner-Skan  equations.  Good  agreement  was  achieved  with 
experimental  data  for  interactions  on  an  adiabatic  wall.  In  the  case  of 
of  an  adiabatic  or  moderately  cooled  wall,  the  enthalpy  profiles  may  be 
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linked  to  the  velocity  profiles,  as  was  done  by  Lees  and  Reeves. 

However,  in  the  case  of  a  iighly  cooled  wall,  the  enthalpy  profile  must 
be  specified  independently  by  use  of  the  energy  equation.  This  was  done 
by  Holden  (Ref.  22),  who  used  a  slightly  modified  Lees  and  Reeves 
(Ref.  21)  approach.  Excellent  agreement  was  obtained  with  experimental 
data  obtained  from  interactions  on  highly  cooled  walls. 

Successful  solutions  to  the  shock-interaction  problem  were  also 
developed  by  Makofski  (Ref.  23)  and  Nielsen  (Ref.  24).  Makofski  obtained 
a  solution  for  laminar  interactions  on  an  adiabatic  wall  by  representing 
the  velocity  profile  by  a  fifth-degree  polynomial  with  two  undetermined 
parameters.  Nielsen  developed  methods  for  laminar  flow  over  both  adia¬ 
batic  and  nonadiabatic  walls. 


EXPERIMENTAL  INVESTIGATIONS 

A  number  of  experimental  investigations  have  been  made  of  shock 
wave-boundary  layer  interactions.  Most  of  the  early  work,  such  as  the 
studies  of  Chapman,  Kuehn,  and  Larson  (Ref.  25),  and  Hakkinen  and  others 
(Ref.  26) ,  was  concerned  with  the  effects  of  the  interaction  on  pressure 
and  skin-friction  distributions  and  boundary-layer  separation.  In  more 
recant  studies,  measurements  of  the  heat-transfer  distributions  through 
interaction  regions  have  been  made.  This  report  correlates  the  results 
of  the  heat-transfer  investigations  in  order  to  obtain  empirical  means 
for  predicting  the  maximum  heat  transfer  and  the  heat-transfer  distri¬ 
bution  in  shock  wave-boundary  layer  interaction  regions.  Results  of  the 
experimental  investigations  described  herein  are  used  in  these  correla¬ 
tions.  The  test  parameters  for  these  investigations  are  summarized  in 
Table  1. 

Levin  and  Fabis'n 


The  thermal  effects  of  shock-wave  impingement  on  a  turbulent 
boundary  layer  were  investigated  in  a  study  conducted  by  Levin  and 
Fabish  (Ref.  27).  Tests  were  conducted  in  a  wind  tunnel  at  free-stream 
Mach  numbers  of  2.95  and  5.02  and  Reynolds  numbers  per  foot  of  7.4  x  10^ 
and  3.6  x  10^,  respectively.  The  boundary  layer  was  generated  by  a 
sharp  leading-edge  flat  plate  instrumented  to  measure  pressure  and 
heating  rates.  Boundary-layer  velocity  profiles  were  measured  upstream 
and  downstream  of  the  interaction.  Turbulent  flow  upstream  of  the 
interaction  region  was  assured  by  means  of  boundary-layer  air  injection 
and  grit.  Air  injection  was  also  used  to  investigate  the  effects  of 
upstream  thickening  of  the  boundary  layer.  The  oblique  shock-wave  gen¬ 
erator  was  varied  in  attitude  from  0  to  12  degrees  and  was  adjusted  so 
that  shock  impingement  occurred  approximately  11.75  inches  from  the 
leading  edge  of  the  flat  plate.  Runs  were  also  made  with  the  shock 
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generator  removed.  Heating  rates  on  the  flat  plate  were  measured  by 
the  use  of  Nichrome  ribbon  heating  elements  that  were  maintained  at  a 
fixed  temperature  above  the  total  temperature  of  the  flow.  The  heating 
rate  was  obtained  from  the  following  relation: 


q  * 


Power  dissipated: 


y  by 


element 


/ 


(Power  conducted 
and  radiated  away 
from  element 


Planform  area  of  element 


) 


Heat-transfer  data  was  presented  in  the  form  of  heat-transfer 
coefficient  obtained  from  the  equation 

^  ^  ^element  ^Awj 

In  some  cases,  shock  impingement  caused  boundary  layer  separation. 
Peak  heat  transfer  always  occurred  within  the  pressure  interaction 
region  and  increased  with  increasing  pressure  rise.  In  the  cases  of 
separated  flow,  peak  heat  transfer  occurred  at  reattachment.  Air 
injection  into  the  boundary  layer  upstream  of  the  interaction  did  not 
significantly  alter  the  heat-transfer  distribution  in  the  interaction 
region. 

Sayano 


In  a  study  similar  to  that  of  Levin  and  Fabish  (Ref.  27),  Sayano 
(Ref.  28)  investigated  the  effects  of  planar-shock  impingement  on  flat 
plates  and  cylinders.  Tests  were  conducted  in  a  20- inch  wind  tunnel. 
Pressure  and  heat-transfer  data  were  obtained  at  free-stream  Mach 
numbers  of  2.41,  3.01,  3.99,  and  5.01.  Free-stream  Reynolds  number 
per  foot  ranged  from  2.56  x  10^  to  4.73  x  10^.  A  boundary-layer  trip 
assured  turbulent  flow  upstream  of  the  interaction.  Runs  were  made 
with  shock  generator  angles  of  8,  15,  and  26  degrees  and  with  the  shock 
generator  removed.  Shock  impingement  occurred  approximately  20  inches 
from  the  leading  edge  of  the  flat  plate  and  the  nose  of  the  cone- cylinder. 

Heat- transfer  measurements  were  made  by  cooling  the  back  surface  of 
the  instrumented  section  and  relating  heat-transfer  coefficient  to  the 
gradient  across  the  wall  of  the  section  by  the  following  expression: 

h  (TR  T0wj  =  b  (T0W  '  Tiw) 
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Radiation  and  conduction  effects  were  neglected,  and  the  accuracy  of 
the  measurements  was  estimated  at  15%.  Pressure  and  heat-transfer 
coefficient  distributions  were  nondimensionalized  by  the  distributions 
obtained  with  the  shock  generator  removed. 

Boundary-layer  separation  occurred  with  shock-generator  angles  of 
15  and  26  degrees.  The  basic  difference  in  shock-impingement  effects 
between  the  flat  plute  and  che  cylindrical  boundary  layers  was  that  the 
interaction  region  was  considerably  larger  and  the  increase  in  pressure 
and  heat-transfer  coefficient  was  greater  on  the  flat  plate.  This 
difference  can  be  attributed  to  crossflow  on  the  cylinder  caused  by  the 
favorable  circumferential  pressure  gradient. 

Peak  heat  transfer  occurred  in  the  same  vicinity  as  peak  pressure 
for  both  the  flat  plate  and  cylindrical  boundary  layers.  In  all  cases, 
the  rise  in  heat-transfer  coefficient  in  the  interaction  region 
increased  with  increasing  shock  strength. 

Magnan  and  Spurlin 

Pianar-shock-impingement  effects  on  laminar  and  turbulent  flat- 
plate  boundary  layers  were  investigated  by  Magnan  and  Spurlin  (Ref.  29). 
Tests  were  conducted  in  a  50-inch-diameter  continuous-flow  wind  tunnel 
at  a  nominal  free-stream  Mach  number  cf  10  and  free-stream  Reynolds 
numbers  per  foot  of  0.5  x  10&  and  2.16  x  10^.  The  impinging  shock  wave 
was  generated  by  a  blunt  leading-edge  flat  plate  set  at  angles  of 
incidence  of  0,  5,  and  10  degrees.  The  position  of  the  shock  generator 
was  adjusted  so  that  shock  impingement  occurred  approximately  22  to 
24  inches  from  the  leading  edge  of  the  flat  plate  in  the  case  of  the 
turbulent  boundary  layer,  and  12  inches  in  the  case  of  the  1  .minar 
boundary  layer. 

The  flat  plate  was  instrumented  to  measure  both  pressure  and 
heating-rate  distributions.  Heating  rates  .'ere  measured  by  the  tran¬ 
sient  thin-skin  method.  The  test  plate,  instrumented  with  thermocouples, 
was  uooled  to  room  temperature  between  runs  before  being  injected  into 
the  air  stream.  The  heating  rate  was  related  to  the  measured  tempera¬ 
ture  response  by  the  expression. 


dT 


W 


dt 


which  neglects  conduction  and  radiation  effects.  Heat-transfer  coef¬ 
ficients  were  obta-,  -.d  from  the  expression 


h  = 


Tw) 
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The  reason  ior  using  total  temperature  rather  than  the  more  conventional 
recovery  temperature  in  this  expression  was  not  given.  If  the  heat- 
transfer  coefficients  of  Ref.  29  are  recalculated  based  on  recovery 
temperature,  the  turbulent  values  are  increased  approximately  20%,  and 
the  laminar  values  approximately  25%.  In  the  correlations,  given  in 
Section  5,  the  peak  heat-transfer  coefficients  (as  obtained  from  Ref.  29) 
are  usually  below  the  general  trend  of  the  data. 

Boundary-layer  separation  occurred  in  the  case  of  shock  impingement 
on  the  laminar  boundary  layer.  Peak  heat-transfer  coefficients  occurred 
in  the  same  region  as  peak  pressures  and  increased  with  increasing  shock 
strength.  Some  difficulty  was  encountered  in  accurately  determining 
pressure  and  heat-transfer  coefficients  from  the  plots  provided  in  the 
report,  particularly  the  values  upstream  of  the  interaction. 

Popinski 

An  experimental  investigation  of  shock  wave-laminar  boundary  layer 
interaction  in  a  two-dimensional  compression  corner  (ramp)  was  conducted 
by  Popinski  (Ref.  30).  Tests  were  performed  in  a  blow-down  wind  tunnel 
at  nominal  free-stream  Mach  numbers  cf  8  and  10  and  Reynolds  numbers 
per  foot  ranging  from  1  x  10^  to  2  x  10^.  For  use  in  the  correlations, 
Mach  numbers  immediately  upstream  of  the  interaction  region  were  cal¬ 
culated  by  assuming  an  isentropic  compression  through  a  leading-edge 
Mach  wave  from  the  free-stream  pressure  to  the  measured  pressure  on  the 
plate. 

Pressure  and  heat-transfer  measurements  were  made  on  a  sharp  leading- 
edge  flat  plate  with  an  adjustable  ramp  that  could  be  set  at  angles 
varying  from  0  to  25  degrees.  The  corner  was  located  13.4  inches  from 
the  leading  edge  of  the  plate. 

Heat-transfer  coefficients  were  determined  by  the  transient  thin- 
skin  method  with  corrections  made  for  spanwise  and  transverse  conduction. 
Radiation  effects  were  neglected.  The  following  relation  was  used: 


The  temperature  response  of  the  flat  plate  was  measured  by  means  of 
thermocouples.  Average  wall  to  total  temperature  ratio  throughout  the 
tests  was  approximately  0.4.  Measured  heat-transfer  coefficients  were 
nondimens ionalized  by  the  theoretical  flat-plate  value  at  the  corner. 
The  distribution  of  pressure  taps  did  not  allow  the  measurement  of 
peak  pressure  on  the  ramp.  For  correlation  purposes,  the  inviscid 
pressure  ratio  calculated  from  the  upstream  Macu  number  and  ramp  angle 
was  used. 
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Boundary-layer  separation  occurred  on  most  of  the  runs.  This  was 
indicated  by  a  plateau  in  the  pressure  distribution  and  a  decrease  in 
heat-transfer  coefficient  below  the  fiat-plate  value  in  the  vicinity  of 
the  corner.  Peak  heat-transfer  coefficients  occurred  on  the  ramp  and 
increased  with  increasing  pressure  rise.  Peak  heat  transfer  was  ob¬ 
served  to  increase  with  decreasing  free-stream  Reynolds  number. 

Needham 


Laminar  interactions  in  a  wedge-type  compression  corner  were  studied 
experimentally  by  Needham  (Ref.  31).  The  object  of  the  investigation  was 
to  determine  if  the  heat- transfer  distribution  could  be  used  as  a  crite¬ 
rion  for  detecting  the  onset  of  separation.  Tests  were  conducted  in  a 
hypersonic  gun  tunnel  at  a  freestream  Mach  number  of  9.7  and  Reynolds 
number  per  foot  of  1.23  x  10^.  The  local  Mach  number  upstream  of  the 
corner  was  determined  by  Shumway  (Ref.  32)  to  be  approximately  8.  Wedge 
angles  were  varied  from  7.6  to  15.3  degrees;  although,  peak  heat-transfer 
measurements  were  presented  only  for  wedge  angles  of  7.6  and  10  degrees. 
The  corner  was  located  6  inches  from  the  leading  edge  of  the  model. 

The  model  was  instrumented  to  measure  pressure  and  heat-transfer 
distributions.  Heating  rates  were  measured  with  t'nin-film  resistance 
thermometer  gages.  Data  were  presented  as  pressure  and  heating  rates 
nondimensionalized  by  flat-plate  values.  Heating  rates  were  converted 
to  heat-transfer  coefficients  by  the  method  developed  by  Shumway. 
Boundary-layer  separation  occurred  at  wedge  angles  of  10  degrees  and 
greater. 

Hastings,  and  Others 

Turbulent  heat  transfer  in  a  wedge  compression-corner  was  inves¬ 
tigated  by  Hastings,  Brown,  and  Atkinson  (Ref.  33).  Tests  were  con¬ 
ducted  in  a  3-  by  4-foot  wind  tunnel  at  a  free-stream  Mach  number  of  4. 
Mach  numbers  upstream  of  the  corner  were  varied  from  2.5  to  4.4  by 
pitching  the  model.  The  corresponding  free-stream  Reynolds  numbers  per 
foot  ranged  from  13.25  x  10^  to  8.33  x  10&. 

The  model,  instrumented  to  measure  pressure  and  heating  rates, 
consisted  of  a  sharp  leading-edge  flat  plate  with  a  15-degree  wedge. 

The  corner  was  located  7.25  inches  from  the  leading  edge  of  the  plate. 
Steady-state  heating  rates,  measured  by  heat  flux  gages,  were  obtained 
by  cooling  the  test  section  of  the  model.  Heat-transfer  data  was 
presented  in  the  form  of  Stanton  number  which  was  obtained  from  the 
expression 
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Stanton  numbers  on  the  ramp  were  referenced  to  conditions  both  upstream 
and  downstream  of  the  corner. 

Boundary- layer  separation  was  not  observed  during  the  testing. 

Although  a  significant  increase  in  Stanton  number  occurred  on  the  ramp, 
heating  rates  both  upstream  and  downstream  of  the  corner  were  lower  than 
had  been  estimated.  The  authors  were  therefore  uncertain  as  to  the 
accuracy  of  the  data.  Heat-transfer  instrumentation  of  the  ramp  did  not 
appear  dense  enough  to  assure  that  the  peak,  heating  rate  would  be  measured. 
Because  of  these  factors,  results  of  this  study  were  used  with  caution  in 
the  present  correlations. 

Becker  and  Korycinski 

Pressure  and  heat- transfer  distributions  in  the  vicinity  of  conical 
flare  compression-corners  were  investigated  in  wind  tunnel  tests  by 
Becker  and  Korycinski  (Ref.  34).  Laminar,  turbulent,  and  transitional 
boundary  layers  were  investigated  at  a  free-stream  Mach  number  of  6.8 
and  Reynolds  numbers  per  foot  from  1.12  x  10?  to  5.6  x  10^.  The  model 
consisted  of  an  ogive-cylinder  with  conical  flare  skirts.  Flare  half¬ 
angles  of  10  and  30  degrees  were  tested.  The  flare-cylinder  junction 
was  15  inches  (10  body  diameters)  from  the  nose  of  the  model. 

Separate  models  were  used  for  pressure  and  heat-transfer  measure¬ 
ments.  Heating  rates  were  determined  by  the  transient  thin-skin  method. 

The  model  was  cooled  to  room  temperature  before  being  injected  into  the 
airstream.  Heat-transfer  data  were  presented  in  terms  of  the  Stanton 
number.  Since  the  Stanton  numbers  were  based  on  free-stream  conditions, 
the  Stanton  number  ratio  was  equal  to  the  heat-transfer-coefficient 
ratio,  which  was  used  in  the  correlations  developed  in  this  report. 

Boundary-layer  separation  occurred  in  the  case  of  laminar  inter¬ 
action  with  a  10-degree  flare.  Heating  rates  in  the  separated  region 
were  approximately  50%  of  the  theoretical  attached  value.  The  measured 
peak-heat  transfer  was  below  the  value  which  was  calculated  by  assuming 
that  a  new  boundary  layer  started  at  the  reattachment  point.  Peak  heat 
transfer  occurred  downstream  of  reattachment. 

In  the  turbulent  case,  a  small  separation  region  existed  with  the 
30-degree  flare,  but  not  with  the  10-degree  flare. 

Extensive  separation  occurred  in  the  case  of  the  transitional 
boundary  layer  with  the  30-degree  flare.  Heating  rates  in  the  separated 
region  were  initially  below  the  attached  value,  but  increased  rapidly 
following  transition.  Peak  heating  rates  occurred  at  reattachment. 
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Schaeffer  and  Ferguson 

Another  experimental  study  of  heat  transfer  in  the  vicinity  of 
conical  flares  was  made  by  Schaeffer  and  Ferguson  (Ref.  35).  Their 
model  differed  slightly  from  that  of  Becker  and  Korycinski  in  having  a 
conical  rather  than  an  ogive  nose.  Tests  were  conducted  in  a  1-foot 
wind  tunnel  at  a  free-stream  Mach  number  of  4.98.  Reynolds  numbers 
per  foot  ranged  from  1.6  x  10&  to  5.4  x  10^  resulting  in  laminar, 
turbulent,  and  transitional  boundary  layers.  Flare  half-angles  of 
10,  17,  24,  and  56  degrees  were  investigated.  The  flare-cylinder 
juncture  was  11.33  inches  from  the  nose  of  the  model,  or  approximately 
6.3  body  diameters. 

The  model,  instrumented  to  measure  pressure  and  heat-transfer 
distributions,  was  cooled  to  120°R  before  being  injected  into  the 
airstream.  Heating  rates  were  determined  by  the  transient  thin-skin 
method  neglecting  conduction  and  radiation  losses.  Accuracy  of  the 
heat-transfer  data  was  estimated  at  ±20%.  Heat-transfer  results  were 
presented  in  the  form  of  free-stream  Stanton  number. 

Separation  occurred  with  the  laminar  and  transitional  boundary 
layers,  but  r.ot  with  the  turbulent  boundary  layer.  The  length  of  the 
separated  region  decreased  with  wall-cooling,  decreasing  flare  angle, 
and  increasing  Reynolds  number.  Minimum  heat  transfer  occurred  slightly 
downstream  of  separation.  Contrary  to  the  results  of  most  of  the  other 
studies  (Ref.  30,  31,  34,  and  36),  heat  transfer  in  the  separated 
laminar  regions  was  not  below  the  attached  value.  Sharp  increases  in 
heat  transfer  occurred  during  transition  and  reattachment.  Peak  heat 
transfer  in  the  pure  laminar  case  occurred  downstream  of  reattachment 
and  was  much  greater  than  the  theoretical  attached  value.  Peak  heating 
in  the  transitional  case  occurred  in  the  reattachment  region.  Peak 
turbulent  heat  transfer  was  predicted  adequately  by  assuming  the  boundary 
layer  started  at  the  flare-cylinder  junction. 

Gulbran,  and  Others 

Three-dimensional  shock  wave-laminar  boundary  layer  interactions 
were  the  subject  of  a  study  made  by  Gulbran  and  others  (Ref.  36).  The 
shock  wave  was  generated  by  a  fin  mounted  perpendicular  to  a  sharp 
leading-edge  flat  plate.  The  test  was  performed  in  a  50-inch-diameter 
wind  tunnel  at  a  free-stream  Mach  number  of  8  and  free--stream  Reynolds 
numbers  per  foot  ranging  from  0.67  x  10&  to  3.0  x  10&.  Parameters 
varied,  in  addition  to  Reynolds  number,  were  fin  leading-edge  sweep 
and  bluntness,  fin  and  plate  angle  of  attack,  and  chordvise  position 
of  the  fin. 
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The  flat  plate,  instrumented  to  measure  pressure  and  heat-transfer 
distributions,  was  cooled  to  O'T  before  being  injected  into  the  airstream. 
The  test  section  of  the  plate  extended  approximately  7  inches  from  thn 
leading  edge.  Heating  rates  were  determined  by  the  transient  thin-skin 
method,  with  corrections  made  for  conduction  losses  when  necessary. 

Data  were  presented  in  the  form  of  pressures  and  heat-transfer  coeffi¬ 
cients  nondimensionalized  by  theoretical  flat-plate  values. 

Results  indicated  that  the  peak  pressure  and  heating  rates  increased 
with  increasing  fin  deflection,  decreasing  sweep  angle,  and  increasing 
leading-edge  diameter.  Peak  heating  rates  were  also  increased  by  moving 
the  fin  forward.  Angle  of  attack  of  the  flat  plate  had  little  effect  on 
the  peak  heating  rates.  Boundary-layer  separation  occurred  in  nearly 
all  test  runs.  Heat-transfer  coefficient  distributions  exhibited  the 
familiar  trend  for  laminar  interactions  of  falling  below  the  attached 
value  in  the  separated  region  and  in  reaching  a  peak  slightly  downstream 
of  reattachment . 

Miller,  and  Others 

As  an  earlier  part  of  the  study  of  Gulbran  (Ref.  36),  Miller  and 
others  (Ref.  37)  investigated  hypersonic  fin-shock  interactions  on  a 
laminar  boundary  layer.  The  unswept,  sharp  leading-edge  fin  was 
mounted  on  a  sharp  leading-edge  flat  plate  instrumented  to  measure 
pressure  and  heat-transfer  distributions.  The  fin  angle  was  varied 
from  -10  degrees  (expansion)  to  +15  degrees  (compression) .  Tests  were 
conducted  in  a  44-inch  wind  tunnel  at  a  free-stream  Mach  number  of  16 
and  Reynolds  numbers  per  foot  ranging  from  0.04  x  10^  to  1.0  x  lO^.  At 
this  high  Mach  number,  interactions  between  the  leading-edge  shock  and 
the  boundary  layer  will  probably  exist;  however,  the  results  of  the 
study  were  included  for  interest  in  the  correlations  presented. 

Heating  rates  were  measured  with  thermistor  heat-flux  gages. 
Heat-transfer  data  were  presented  in  the  form  of  free-stream  Stanton 
number.  Measured  pressures  were  nondimensionalized  by  the  free-stream 
static  pressure. 

Heat-transfer  results  indicated  that  boundary- layer  transition 
occurred  at  compression  angles  greater  than  10  degrees.  Pressure  and 
heat-transfer  distributions  did  not  indicate  the  presence  of  separation. 

Experimental  investigations  of  heat  transfer  on  conical  flares  and 
deflected  flaps  also  were  made  by  utilizing  free-flight  testing  (Ref.  38 
and  39).  These  results  were  not  used  in  the  present  correlations 
because  of  the  .‘elative  inaccuracy  of  measurements  made  during  free 
flight  as  compared  with  wind-tunnel  testing. 
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EXISTING  EMPIRICAL  RELATIONS  FOR  PEAK 
HEAT  TRANSFER 

Several  of  the  experimental  investigators  used  their  test  results 
to  develop  empirical  expressions  for  determining  the  peak  heat  transfer 
in  shock  interaction  regions.  Most  of  these  expressions  were  developed 
with  a  limited  amount  of  data  for  a  specific  configuration.  Relations 
that  were  presented  in  some  of  the  studies  discussed  in  the  previous 
section  are  described  below. 

Levin  and  Fabish 

An  expression  for  peak  Stanton  number  in  the  interaction  region  of 
an  impinging  planar  shock  wave  on  a  turbulent  flat-plate  boundary  layer 
was  developed  by  Levin  and  Fabish,  from  their  experimental  data  (Ref.  27). 
They  suggested  the  following  relation  between  the  peak  Stanton  number 
and  the  peak  pressure  ratio: 


The  constant  K^  has  a  value  from  approximately  1.4  to  1.8;  K2  is 
approximately  1.45;  and  the  reference  value  (P3/Pl)g  is  a  function  of 
the  mass  flow  ratio  (pV)3/(pV)i. 


Sayano 


A  relation  between  peak  heat  transfer  and  peak  pressure  was  pre¬ 
sented  in  Sayano's  study  (Ref.  28).  The  expression  for  turbulent 
interactions  arising  from  planar  shock  impingement  on  flat  plates  or 
cylinders  is 
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The  peak  pressure  ratio  was  also  related  to  the  incident  shoct-: 
strength;  the  resulting  expressions  for  peak  heat  transfer  are 


Flat  Plate  Impingement:  —  = 

nFP 


Cylinder  Impingement:  r —  = 

FP 


Popinski 

Popinski  performed  peak  heat-transfer  correlations  using  both  his 
laminar  data  and  the  results  of  other  investigations.  The  following 
expression  was  obtained  for  the  case  of  laminar  interactions  in  two- 
dimensional  compression  corners r 


For  separated  turbulent  interactions  caused  by  shock  impingement,  he 
suggested  the  following  relation  between  peak  heat  transfer  and  plateau 
pressure 


For  the  case  of  turbulent  interactions  without  separation,  peak  heat 
transfer  was  related  to  pressure  ratio  and  Reynolds  number  as  follows: 


Both  of  the  turbulent  relations  were  obtained  from  the  data  of  Levin 
and  Fabish  (Ref.  27)  and  Sayano  (Ref.  28). 
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In  addition  to  peak  heat-transfer  correlations,  Popinski  developed 
empirical  expressions  for  several  other  interaction  parameters  such  as 
plateau  pressure,  interaction  length,  and  separation  length. 

Gulbran,  and  Others 


The  results  of  Gulbran  and  others  (Ref.  36)  for  laminar  fin-shock 
interactions  indicated  that  the  governing  variable  for  the  peak  heat 
transfer  was  che  rise  from  plateau  pressure  to  peak  pressure.  They 
proposed  the  following  empirical  expression: 


In  this  particular  correlation,  PPLj^p  is  f°r  case  of  a  sharp 

leading-edge  unswept  fin  at  7.5  degrees  deflection.  The  plateau 
pressure,  P?l,  is  calculated  by  the  method  of  Erdos  and  Pallone  (Ref.  40), 
and  the  final  pressure,  Pp,  is  obtained  from  oblique  shock  theory. 

Shuraway 

Results  of  several  experimental  studies  were  used  by  Shumway 
(Ref.  32)  to  obtain  a  relation  for  the  peak  heat  transfer  in  shock 
wave  boundary  layer  interaction  regions.  He  suggested  the  following 
linear  relation  between  peak  heat-transfer  coefficient  and  peak  pressure 
in  the  interaction  region: 


h 

u 


P 

u 


This  expression  was  obtained  by  using  data  from  laminar  and  turbulent 
interactions  generated  by  compression  corners,  impinging  shock  waves, 
and  fin  shocks.  The  subscripts  d,  downstream,  and  u,  upstream,  used 
by  Shumway  are  equivalent  to  the  subscripts  PK,  peak,  and  FP,  flat  plate, 
which  are  used  in  the  studies  mentioned  previously. 

A  similar  expression  can  be  obtained  by  order-of-magnitude  analysis, 
as  shown  by  Giles  and  Thomas  (Ref.  41).  Because  of  certain  assumptions 
made  in  the  analysis,  the  authors  recommend  that  this  result  be  used 
only  in  a  qualitative  sense. 
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For  cases  in  which  shock  interaction  causes  transition,  Shumway 
suggests  a  relation  based  on  the  following  laminar  and  turbulent 
expressions  for  Stanton  number: 


St  0.332  Pr“2/3  Re”°‘5 

LAM 


-2/*}  -0  2 

StTURB  =  °*°283  Pr  Re 


The  ratio  of  these  expressions,  taken  across  a  shock  wave,  gives 
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u 


(PV) 


0.8 


0.089 


(PV) 
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Section  4 

PRESSURE  AND  HEAT  TRANSFER  DISTRIBUTIONS 


Schematic  representations  indicating  the  nomenclature  used  to 
define  the  geometry  of  the  impinging-shock,  compression-corner,  and 
fin-shock  interactions  are  given  in  Fig.  10.  Qualitative  distributions 
of  pressure  and  heat-transfer  coefficient  in  interaction  regions  are 
shown  in  Fig.  11  and  12,  respectively.  The  reference  position,  shown 
in  these  figures,  represents  the  point  of  geometric  shock  impingement 
in  the  case  of  the  impinging  shock;  the  corner  for  wedge  or  conical 
flare  interactions;  and  in  the  case  of  the  fin  shock,  the;  location 
where  the  shock  is  incident  on  the  surface. 

When  analyzing  shock  wave-boundary  layer  interactions,  it  is  first 
necessary  to  know  whether  the  interaction  is  laminar  (transition 
downstream  of  interaction);  turbulent  (transition  upstream  of  inter¬ 
action);  or  transitional  (transition  within  interaction  region). 
Unfortunately,  no  completely  reliable  way  is  known  for  predicting 
whether  or  not  the  latter  will  occur. 

It  is  also  necessary  co  determine  If  the  boundary  layer  is  separated 
or  attached.  Some  relations  have  been  developed  for  predicting  the 
occurrence  of  separation.  Incipient  separation  is  defined  as  the  con¬ 
dition  under  which  flow  first  begins  to  separate.  Separation  will  occur, 
therefore,  if  the  overall  pressure  rise  associated  with  the  interaction 
is  greater  than  the  pressure  rise  required  for  incipient  separation. 

This  pressure  is  sometimes  termed  the  critical  pressure. 

Using  the  results  of  a  study  of  shock  impingement  on  laminar 
boundary  layers,  Hakkinen  and  others  (Ref.  26)  obtained  the  following 
relation  between  the  pressure  coefficient  required  for  incipient  sep¬ 
aration,  and  the  pressure  coefficient  at  the  separation  point, 

CpSEP* 


The  following  empirical  expression  for  separation  pressure  in 
laminar  interactions  was  developed  by  Popinski  (Ref.  30)  and  is  based 
on  tests  of  several  configurations: 


C  =  0.91  (M^  - 
SEP  '  1 


l\-°'287  R,-174 

/  °° 
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(a)  Shock  Impingement  on  Flat 
Plate  or  Cylinder. 


(b)  Compression  Corner  With  Attached 
Boundary  Layer. 


(c)  Compression  Corner  With  Separated 
Boundary  Layer. 


FIG.  10.  Representation  of  Nomenclature  Used 
to  Describe  Interaction  Parameters. 
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Combining  the  two  equations  gives 


C 


P 


INC 


1.82 


-0.287 


4 


x 


0 


The  Reynolds  number  in  this  equation  is  calculated  at  the  beginning 
of  the  interaction  region,  xq.  However,  if  xo  is  large  compared  to  the 
length  of  interaction  region,  the  Reynolds  number  may  be  evaluated  at 
xref  without  having  an  appreciable  effect  on  Cp^^,.  Although  the  rela¬ 
tion  between  Cpt„_  and  Cp(._0  is  based  on  data  from  impinging  shock 
INC  Sfcr 

interactions,  if  should  be  valid  for  all  configurations  because  of  the 
free-interaction  concept,  which  states  that  the  upstream  portion  of  the 
interaction  is  independent  of  the  means  of  generating  the  pressure 
gradient  (Ref.  25). 

The  following  expression  for  incipient  separation  in  turbulent 
boundary  layers  was  developed  by  Popinski  who  used  compression  corner 
data: 


r 
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As  in  t:;e  case  of  laminar  separation,  this  expression  should  be  appli¬ 
cable  to  other  configurations  because  of  the  concept  of  free  interaction. 


PRESSURE  DISTRIBUTION 

Figure  II  indicates  the  importance  of  the  state  of  the  boundary 
layer  on  the  pressure  distribution  through  an  interaction  region.  If 
the  boundary  layer  remains  attached,  the  pressure  rises  smoothly  from 
the  undisturbed  value  to  the  peak  value.  The  pressure  rise  in  a 
laminar  interaction  will  take  place  over  a  considerably  longer  distance 
than  it  will  in  the  case  of  a  turbulent  interaction.  This  is  because 
of  the  lower  kinetic  energy  of  the  laminar  boundary  layer  allowing  the 
pressure  disturbance  to  propagate  farther  upstream.  If  the  boundary 
layer  separates,  the  pressure  increases  at  the  separation  and  reattach¬ 
ment  points,  and  a  pressure  plateau  forms  in  the  separated  region.  In 
separated  turbulent  flow,  the  plateau  will  often  appear  3S  an  inflection 
point  because  the  separated  region  is  of  relatively  short  length.  In 
addition  to  the  separation  and  reattachment  pressure  rises,  pressure 
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rises  have  been  observed  in  transition  regions  of  separated  flows. 

A  detailed  discussion  of  the  effects  of  the  state  of  the  boundary  layer 
on  the  pressure  distribution  may  be  found  in  the  work  of  Chapman,  Kuehn, 
and  Larson  (Ref.  25).  Plateau  pressure,  peak  pressure,  and  the  inviscid 
final  pressure  are  used  in  the  present  heat-transfer  correlations. 

Plateau  Pressure 


A  region  of  constant  pressure,  designated  the  plateau  pressure, 
has  been  observed  to  exist  in  separated  regions.  Plateau  pressure  is 
governed  by  the  concept  of  free  interaction  and  is  a  function  only  of 
the  upstream  Mach  number  and  Reynolds  number.  Analysis  (Ref.  25)  has 
shown  that 


1/2 


With  the  use  of  the  following  equations  for  laminar  and  turbulent 
local  skin  friction  (Ref.  42)  , 


=  0.664  Re 


-1/2 


x 


=  0.0576  Re 


-1/5 


TURB 


expressions  for  laminar  and  turbulent  plateau  pressures  are 


TURB 
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Several  values  of  Che  constants  Kl  and  Kj  were  proposed  from  the 
results  of  various  studies  (Ref.  26,  40,  and  43).  For  the  present 
correlations,  the  values  Kl  =  1.6  and  K?  =  1.7  are  used,  as  given  in 
the  study  of  Wuerer  and  Clayton  (Ref.  43).  Plateau  pressures  calculated 
by  using  the  preceding  equations,  and  experimental  plateau  pressures 
obtained  from  the  studies  used  in  the  present  correlations  are  compared 
in  Fig.  13.  The  agreement  is  satisfactory,  with  the  exception  of  two 
data  points  (Ref.  35).  Some  of  the  scatter  can  be  attributed  to  the 
difficulty  of  identifying  the  plateau  pressure  in  several  of  the  exper¬ 
imental  pressure  distributions. 


EXPERIMENTAL  PLATEAU  PRESSURE ,  P^tjy  /  P„ 

FIG.  13.  Comparison  of  Theoretical  and  Experimental 
Plateau  Pressures. 


In  the  case  of  certain  configurations,  the  plateau  pressure  may  be 
used  to  estimate  the  length  of  the  separated  region.  The  flow- 
deflection  angle  in  the  separated  region,  0',  can  be  obtained  from 
the  plateau  pressure  and  upstream  Mach  number  by  using  inviscid  shock 
relations.  If  either  the  separation  or  the  reattachment  points  are 
known  on  a  compression  corner,  the  separation  length  may  be  determined 
geometrically  by  using  the  flow-deflection  angle. 
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Peak  and  Final  Pressures 


After  passing  through  the  shock  wave,  the  pressure  distribution 
reaches  a  peak  value  which  may  or  may  not  exceed  the  theoretical  inviscid 
pressure-rise  to  the  final  pressure  behind  the  interaction.  Analytical 
means  for  predicting  the  peak  pressure  associated  with  shock  wave¬ 
boundary  layer  interactions  have  been  developed  (Ref.  44  and  45).  In 
most  cases,  however,  the  theoretical  final  pressure  is  a  good  approxi¬ 
mation  of  the  peak  pressure.  In  Fig.  14,  the  corresponding  final 
pressures  are  plotted  against  the  experimental  peak  pressures  obtained 
from  the  studies  used  in  the  correlations  presented.  Final  pressures 
are  calculated  by  using  the  applicable  inviscid  shock  theories  as 
described  below.  A  more  detailed  description,  including  the  equations 
involved,  is  given  in  the  Appendix. 


FIG.  14.  Comparison  of  Theoretical  Final 
Pressure  With  Experimental  Peak  Pressure. 
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Shock  Impingement  on  Flat  Plate  or  Cylinder.  For  the  case  of  shock 
impingement  and  reflection  (Fig.  10(a)),  the  pressure  rise  across  the 
incident  shock  is  calculated  from  the  upstream  Mach  number  and  the  shock- 
generator  deflection  angle  with  the  use  of  the  oblique  shock  relations 
(Ref.  46).  The  pressure  rise  across  the  reflected  shock  to  the  final 
pressure  is  calculated  from  the  Mach  number  behind  the  incident  shock 
with  a  turning  angle  equal  to  the  shock-generator  angle  so  that  the  flow 
is  turned  back  parallel  to  the  surface.  The  reflected  shock  is  assumed 
to  be  oblique  for  the  flat-plate  interaction  and  conical  for  the  cylinder 
interaction. 

Wedge  Compression-Corner.  The  pressure-rise  across  an  attached- 
wedge  interaction  (Fig.  10(b))  is  calculated  from  the  upstream  Mach 
number  and  the  wedge  angle.  If  the  flow  separates,  the  pressure  rise 
is  assumed  to  take  place  through  a  double  compression.  The  initial 
pressure  rise  is  assumed  to  be  the  plateau  pressure  calculated  from  the 
plateau  pressure  equations  given  on  page  28.  The  flow-deflection  angle 
is  calculated  from  the  plateau  pressure  and  the  upstream  Mach  number. 

The  rise  from  plateau  pressure  to  final  pressure  is  then  calculated  by 
using  the  Mach  number  behind  the  separation  shock  and  by  assuming  that 
the  flow  turns  parallel  to  the  wedge  surface.  Since  the  plateau  pres¬ 
sure  is  based  on  skin  friction,  this  is  not  truly  an  inviscid  pressure 
rise.  The  double-compression  pressure  rise  may  be  approximated  by 
assuming  that  the  separation  angle  is  half  the  size  of  the  wedge  angle. 
Figure  i5  shows  how  the  final  pressure  varies  with  separation  angle  for 
a  20-degree  wedge  at  Mach  4.  Since  the  separation  angle  is  usually 
small  compared  to  the  overall  turning  angle,  this  approximation  will 
yield  conservative  results. 

Conical  Flare  Compression-Corner.  The  pressure  rise  across  attached 
or  separated  conical-flare  interactions  was  calculated  in  the  same  man¬ 
ner  as  that  of  the  wedge  interactions,  with  the  exception  that  empirical, 
conical  shock  relations  (Ref.  47)  were  used  rather  than  the  oblique  shock 
relations . 

Fin  Shock.  The  pressure  rise  across  a  fin-generated  shock  wave 
(Fig.  10(d))  was  calculated  from  the  upstream  Mach  number  and  the  fin- 
deflection  angle  by  using  the  oblique  shock  relations. 

In  some  cases,  detached  shock  waves  were  predicted,  and  as  a  result 
solutions  for  final  pressures  were  not  obtained. 
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9  e  20  DEG 


FIG.  15.  Variation  of  Inviscid  Pressure 
Ratio  With  Separation  Angle. 


HEAT-TRANSFER  COEFFICIENT  DISTRIBUTION 

As  would  be  expected,  the  heat-transfer  coefficient  distribution, 
shown  qualitatively  in  Fig.  12,  is  strongly  dependent  on  the  state  of 
the  boundary  layer  in  the  interaction  region  as  is  the  pressure  dis¬ 
tribution.  It  the  boundary  layer  remains  attached  through  the  inter¬ 
action,  the  heat-transfer  coefficient  distribution  follows  the  same 
trend  as  that  of  the  pressure  distribution,  increasing  relatively 
smoothly  from  the  beginning  of  the  interaction  to  the  peak  value.  If 
the  boundary  layer  separates,  or  if  transition  occurs  in  the  interaction 
region,  the  heat-transfer  coefficient  distribution  becomes  considerably 
more  involved.  The  areas  of  interest  in  the  heat-transfer  coefficient 
distributions  are  heat  transfer  in  the  separated  region  and  peak  heat 
transfer. 
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Heat-Transfer  Coefficients  in  Separated  Regions 

Analysis  of  separated  flow  is  complicated  because  it  is  difficult 
to  obtain  repeatable  experimental  results  in  consecutive  runs,  let  alone 
to  compare  results  from  different  facilities  (Ref.  41).  Some  of  the 
trends  which  have  been  noticed  are  discussed  below. 

Laminar  Separation.  There  is  a  definite  decrease  in  heat-transfer 
coefficient  in  separated  laminar  regions.  Chapman  (Ref.  48)  performed 
a  theoretical  analysis  which  indicated  that  the  average  laminar  separated 
heating  rate  was  56%  of  the  attached  laminar  value.  This  is  in  agreement 
with  several  experimental  studies  (Ref.  30,  34,  49,  50,  and  51)  that 
show  the  separated  value  to  be  approximately  half  the  attached  value. 

Turbulent  Separation.  No  reliable  means  for  predicting  heat  trans¬ 
fer  in  separated  turbulent  flows  is  known.  Chapman  (Ref.  48)  conducted 
a  limited  analysis  showing  heating  rates  above  the  attached  values, 
with  the  rate  of  heating-increase  decreasing  as  the  Mach  numbers 
increased.  Most  experimental  studies  show  the  separated  turbulent  heat- 
transfer  to  be  above  the  attached  value  (Ref.  28,  51,  and  52);  however, 
heating  rates  below  the  attached  value  have  been  measured  (Ref.  49). 
Semi-empirical  relationships  for  the  calculation  of  Stanton  number  in 
regions  of  turbulent  separated  flows  were  developed  by  Spalding 
(Ref.  53)  that  satisfactorily  correlate  with  experimental  results. 

Peak  Heat-Transfer  Coefficient 


Peak  heat  transfer  in  shock  wave-boundary  layer  interaction  regions 
is  usually  of  the  utmost  interest.  Heat-transfer  coefficient  increases 
of  10  tj  '•  the  undisturbed  value  (the  value  which  would  exist  with  no 
interaction)  are  common;  in  cases  where  the  interaction  causes  transi¬ 
tion,  increases  of  over  30  times  the  undisturbed  value  have  been 
measured.  The  seriousness  of  a  high  local  heat-transfer  coefficient 
is  related  to  the  thermal  properties  of  the  material  involved.  If 
the  surface  material  has  a  high  thermal  conductivity,  heat  con¬ 
duction  away  from  the  region  of  high  heating  rate  will  provide  relief 
from  thermal  gradients.  If  the  thermal  conductivity  is  low,  however, 
heat  conduction  will  be  much  less,  and  a  hot  spot  will  develop. 

Peak  measured  heating  rates  have  usually  been  observed  in  the 
vicinity  of  reattachment  where  the  boundary  layer  is  thinnest,  or  if 
-.paration  does  not  occur,  immediately  downstream  of  shock  interaction. 
Increased  mixing  rates  in  the  interaction  region  contribute  to  the 
increase  in  heat  transfer.  Some  success  has  been  achieved  in  predicting 
peak  heat  transfer  by  assuming  that  a  new  boundary  layer  starts  at  shock 
impingement  or  at  a  compression  corner  for  attached  flows,  and  at  re¬ 
attachment  for  separated  flows.  As  mentioned  previously,  several  cor¬ 
relations  of  peak  heat  transfer  have  been  made  for  specific  configurations. 
In  Section  5,  experimental  peak  heat-transfer  coefficients  for  several 
configurations  are  correlated  with  various  parameters. 
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Little  work  has  been  done  on  the  decay  of  heat-transfer  coefficient 
distribution  downstream  of  shock  wave-boundary  layer  interactions. 


INTERACTION  GEOMETRY 
Interaction  Lengths 

Several  correlations  of  the  interaction  lengths  that  are  required 
to  define  the  pressure  and  heat-transfer  distributions  were  developed 
by  Popinski  (Ref.  30).  Expressions  are  presented  for  interaction 
lengths  upstream  and  downstream  of  both  laminar  and  turbulent  boundary- 
layers.  The  correlations  are  d5vided  into  impinging-shock  data  and 
compression-corner  data,  and  in  most  cases  they  represent  extremely 
limited  Mach  number  ranges.  The  correlating  parameters  are  Mach  number, 
Reynolds  number,  and  plateau  pressure.  Further  data  of  this  nature  are 
required  before  expressions  for  general  application  can  be  obtained. 

Separation  Length 

Estimates  of  length  of  the  separated  region  for  laminar  (Ref.  41) 
and  for  turbulent  (Ref.  30)  boundary  layers  have  been  made  by  combining 
impinging-shock  and  compression-corner  data.  Although  the  separation 
length  is  defined  as  the  distance  from  separation  to  reattachment,  the 
turbulent  correlation  is  for  the  distance  from  separation  to  the  ref¬ 
erence  line.  Turbulent  separation  is  usually  small  enough  so  that  the 
pressure  and  heat-transfer  coefficient  distributions  can  be  approxi¬ 
mated  by  linear  increases  over  the  interaction  region.  The  expressions 
for  laminar  and  for  turbulent  separation  lengths  are  shown  below. 
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RECOVERY  FACTOR 

Another  important:  aerodynamic-heating  parameter  which  is  affected 
by  shock  wave-boundary  layer  interaction  is  the  recovery  factor,  defined 
as 


r 


The  adiabatic  wall  or  recovery  temperature  is  obtained  from  the  free- 
stream  temperature  and  Mach  number  by  the  relation 

1  +  (r  ^  M») 

The  recovery  factor  is  thus  related  to  the  aerodynamic-heating  rate, 
as  shown  by  the  expression 


’  •  h  (taw  -  Tw) 


Experimental  studies  have  shown  that  the  recovery  factor  In  inter¬ 
action  regions  is  higher  than  the  normal  value,  which  is  approximately 
0.85  for  laminar  flow  and  0.89  for  turbulent  flow.  Means  for  predicting 
this  increase  do  not  exist.  Since  most  of  the  studies  determined  the 
adiabatic  wall  temperature  experimentally,  effects  of  the  variable 
recovery  factor  were  automatically  included  in  the  data  analysis.  Some 
of  the  studies,  however,  assumed  the  appropriate  recovery  factor  and 
calculated  the  adiabatic  wall  temperature.  As  a  result,  a  portion  of 
the  appaient  increase  in  heat-transfer  coefficient  was  actually  an 
increase  in  recovery  factor.  The  errors  caused  by  this  approach  more 
than  likely  were  small  compared  to  the  overall  experimental  accuracy. 
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Section  5 

PEAK  HEAT-TRANSFER  CORRELATIONS 


Examination  of  existing  experimental  data  led  to  the  conclusion 
that  the  only  thermal  quantity  in  the  interaction  region  for  which 
meaningful  correlations  can  be  made,  at  the  present  time,  is  peak  heat: 
transfer.  From  an  aerodynamic  heating  standpoint,  this  is  probably  the 
most  important  quantity  in  the  interaction  region.  However,  because  of 
the  extreme  variations  in  heating  rate  that  occur  through  the  inter¬ 
action,  a  complete  analysis  of  a  specific  problem  would  require  some 
knowledge  of  the  heating  rate  distribution. 

For  this  study,  correlations  of  experimental  peak  heat-transfer 
coefficients  were  made  with  the  following  parameters: 

1.  Peak  pressure  measured  in  the  interaction  region 

2.  Theoretical  final  pressure 

3.  Theoretical  plateau  pressure 

4.  Reynolds  number  based  on  distance 

5.  Reynolds  number  based  on  boundary- layer  thickness 

6.  Mach  number 

Data  were  utilized  from  tests  performed  on  the  following  configura¬ 
tions  : 

1.  Planar  shock  impingement  on  flat  plate  (Ref.  27,  28,  and  29) 

2.  Planar  shock  impingement  on  cylinder  (Ref.  28) 

3.  Compression  corner,  wedge  (Ref.  30,  31,  and  33) 

4.  Compression  corner,  conical  flare  (Ref.  34  and  35) 

5.  Fin  mounted  on  flat  plate  (Ref.  36  and  37) 

These  data  represent  free-stream  Mach  numbers,  ranging  from  2.4  to  16 
and  Reynolds  numbers  based  on  distance  from  7.9  x  10^  to  8.09  x  10&. 
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Pressures  and  heat-transfer  coefficients  are  referenced  to  flat- 
plate  values,  which  were  obtained  either  experimentally  or  analytically, 
depending  on  how  the  original  data  were  presented.  The  flat-plate  values 
represent  the  conditions  that  would  exist  at  the  reference  line  (Fig.  10) 
in  undisturbed  flow.  When  practical,  empirical  expressions  were  obtained 
by  fitting  a  straight  line  through  the  data  points  using  the  method  of 
least  squares.  Parameters  used  in  the  correlations  are  summarized  in 
Table  2.  Results  of  the  correlations  follow. 


PEAK  PRESSURE 

The  first  peak  heat-transfer  correlations  attempted  in  this  study 
were  with  the  experimental  peak  pressures.  Results  of  these  correla¬ 
tions  indicated  that  the  relation  between  peak  heat-transfer  coefficient 
and  peak  pressure  is  independent  of  configuration,  at  least  for  the  five 
configurations  considered.  The  results  also  indicated  that  the  correla¬ 
tions  can  be  separately  categorized  as  to  whether  the  boundary  layer  is 
completely  laminar,  completely  turbulent,  or  undergoes  transition  in 
the  interaction  region.  Plots  of  peak  heat-transfer  coefficient  versus 
peak  pressure  are  shown  for  the  laminar,  turbulent,  and  transitional 
interactions  in  Fig.  16,  17,  and  18,  respectively.  Considerable  scatter 
exists  in  the  laminar  data,  but  the  turbulent  data  correlates  fairly 
well.  Only  a  limited  amount  of  transitional  data  are  available  for 
correlation  purposes. 

Empirical  expression  describing  the  correlations  are 


Laminar: 


Turbulent : 


Transition : 


(based  on 
limited  data 
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TABLE  2.  Summary  of  Test  Results  and  Correlation  Parameters 


Ref. 

Conf ig- 
uration 

M. 

1 

Experimental  1 

Calculated 

Re*nEF 
x  10~6 

r*6 

x  10~4 

Boundary 

layer 

D 

m 

!pk 

PFP 

m 

1 

27 

Flat 

2.95 

6 

2.98 

2.  78 

2.39 

n 

2.34 

(a) 

7.79 

11.9 

Turbulent 

pl3te 

2.95 

6 

2.98 

2.44 

2.35 

2.34 

(a) 

7.79 

.  .  >  . 

Turbulent 

2.95 

9.5 

2.98 

3.19 

3.46 

3.59 

(a) 

8.03 

11.6 

Turbulent 

2.95 

9.5 

2.98 

3.16 

3.56 

3.59 

(a) 

.... 

Turbulent 

2.95 

12 

2.98 

4.11 

5.50 

4.  75 

8.09 

11.3 

TurK  lent 

2.95 

12 

2.98 

3  79 

5.25 

4.75 

BE1 

8.09 

.... 

T-  rbulent 

2.95 

12 

2.98 

3.60 

5.25 

4.75 

2.31 

8.09 

.... 

Turbulent 

5.02 

7.4 

4.77 

4.70 

4.15 

4.60 

(a) 

4.31 

.... 

Turbulent 

5.02 

9.7 

4. 77 

5.63 

5.86 

6.  75 

(a) 

4.24 

.  .  •  • 

Turbulent 

5.02 

9.7 

4.77 

5.35 

6.16 

6.75 

(a) 

4.24 

.... 

Turbulent 

5.02 

11 

4.77 

6.05 

7.22 

8.24 

(a) 

4.24 

.... 

Turbulent 

5.02 

11 

4.77 

6.06 

7.25 

8.24 

(a) 

4.24 

.... 

Turbulent 

5.02 

11 

4.77 

5.67 

6.81 

8.24 

(a) 

4.31 

.... 

Turbulent 

5.02 

12 

4.  77 

7  39 

8.41 

9.52 

4.24 

.... 

Turbulent 

5.02 

12 

4.77 

7.24 

8.47 

9.52 

m 

4.24 

.... 

Turbulent 

CO 

Flat 

2.41 

8 

2.41 

2.0 

2.4 

2.55 

(a) 

7.9 

12.4 

Turbulent 

plate 

2.41 

15 

2.41 

3.4 

4.3 

5.1b 

1.95 

7.9 

12.4 

Turbulent 

3.01 

8 

3.01 

2.2 

3.0 

3.03 

(a) 

7.55 

12.0 

Turbulent 

3.01 

IS 

3.0. 

4.6 

6.2 

6.58 

2.31 

7.55 

12.0 

Turbulent 

3.01 

26 

3.fi 

6.0 

J.  5 

(b) 

2.31 

7.55 

12.0 

Turbulent 

3.99 

8 

3.99 

3.2 

4.0 

4.06 

(a) 

7.05 

11.3 

Turbulent 

'.99 

15 

3.99 

6.8 

11.0 

10.19 

2.99 

7.05 

11.3 

Turbulent 

5.01 

15 

3.01 

. 

9.2 

15.0 

15.53 

3.86 

5.35 

Turbulent 

28 

"ylinder 

2.41 

8 

2.41 

1.6 

2.1 

2.0 

(a) 

7.77 

Turbulent 

2.41 

15 

2.41 

2.2 

3.0 

3.5 

1.95 

7.77 

Turbulent 

3.01 

8 

3.01 

2.2 

2.6 

2.3 

(a) 

7.17 

11.5 

Turbulent 

3.01 

15 

3.01 

3.5 

5.4 

4.7 

2.32 

7.17 

11.5 

Turbulent 

3.99 

8 

3.99 

2.6 

2.7 

3.0 

(a) 

6.97 

11.2 

Turbulent 

3.99 

15 

3.99 

4.8 

■SI 

2.99 

6.97 

11.2 

Turbulent 

5.01 

8 

5.01 

3.3 

ksh 

■ 

(a) 

3.88 

7.0 

Turbulent 

5.01 

15 

5.01 

7.0 

10.7 

11.4 

3.96 

3.88 

7.0 

Turbulent 

29 

rmm 

m 

0 

C 

5 

l 

4.32 

7.67 

Turbulent 

ISM 

ESS 

5 

I  10.18 

7 

10 

7.95 

5.04 

8.66 

Turbulent 

10 

13 

26 

30.35 

5.22 

8.87 

Turh  jlent 

K£J 

10 

Isa 

20.5 

26 

30.35 

0.  79 

0.45 

Laminar 

iO 

Wedge 

8.13 

9.5 

7.8 

9.4 

5.17 

1.67 

2.30 

0.76 

Laminar 

8.12 

15.0 

7.8 

19.8 

10.52 

1.68 

2.11 

0.73 

Laminar 

8.07 

13.9 

7.8 

17.6 

9.18 

1.73 

1.59 

0.63 

Laminar 

8.02 

9.4 

7.8 

17.1 

5.02 

1.  79 

1.11 

0.53 

Laminar 

8.12 

18.  5 

7.4 

30 

15.23 

1.69 

0.72 

Laminar 

9.7 

wm 

9.1 

10.8 

7.32 

1.94 

1.71 

Laminar 

9. 7 

13.  7 

9.7 

22.1 

12.54 

1.94 

1 .  ‘  8 

0.65 

Laminar 

9.68 

10.0 

9,2 

16.  j 

7.30 

1.97 

1.-5 

0.61 

Laminar 

0.7 

10.0 

5.5 

23.6 

20.93 

195 

i.63 

0.64 

Laminar 

I 
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TABLE  2.  (Continued) 


Ref. 

Config¬ 

uration 

M 

'no 

e, 

deg 

M1 

Experimental 

Calculated 

Re*REF 
x  10-6 

r*6 

x  10-* 

Boundary 

layer 

_PK 

hFP 

PPK 

PFP 

PF 

PFP 

p 

PL 

PFP 

31 

W->dge 

B9 

7.6 

8 

4.0 

3.9 

3.8 

(a) 

KE9 

0.47 

Laminar 

Hi 

10.0 

8 

5.9 

5.15 

5.45 

1.83 

Hi 

0.47 

Laminar 

33 

Wedge 

2.495 

Bfl 

2.495 

2.5 

2.5 

2. 4£ 

(a) 

8.0 

rh 

Turbulent 

3.23 

3.23 

3.0 

3.05 

3.00 

(a) 

7.98 

Turbulent 

3.59 

tsm 

3.59 

3.13 

3.35 

3.31 

(a) 

7.15 

11.44 

Turbulent 

3.97 

15 

3.97 

2.9 

3.7 

3.67 

(a) 

5.99 

9.92 

Turbulent 

4.36 

15 

4.  36 

3.45 

4.1 

4.06 

(a) 

5.03 

8.47 

Turbulent 

34 

Conical 

6.8 

6.8 

3.2 

2.3 

2.37 

1.54 

1.9 

m 

Laminar 

flare 

6.8 

6.8 

3.3 

3.2 

2.37 

(a) 

7.0 

Hi 

Turbulent 

6.8 

IB 

6.8 

13.3 

25.6 

19.49 

(a) 

6.3 

isa 

Turbulent 

6.8 

H 

6.8 

43.8 

24.9 

19.49 

1.58 

1.4 

HI 

Transition 

35 

Conical 

4.98 

10 

4.98 

4.3 

2.56 

1.  75 

1.36 

na 

0.62 

Laminar 

flare 

4.98 

10 

4.98 

11.0 

2.74 

1.75 

1.26 

13 

Transition 

4.98 

24 

4.98 

29.0 

8.64 

6.0 

1.26 

Transition 

4.98 

24 

4.98 

5.9 

7.94 

6.0 

3.85 

5.15 

8.82 

Turbulent 

4.98 

56 

4.98 

9.2 

9.5 

(b) 

2.85 

5.15 

8.82 

Turbulent 

36 

Fin 

8.0 

7.5 

3.0 

2.9 

Laminar 

8.0 

7.5 

.  .  . 

5.4 

2.8 

Laminar 

8.0 

7.5 

.  .  . 

6.9 

2.8 

Laminar 

8.0 

15.0 

.  .  . 

6.2 

6 

Laminar 

8.0 

15.0 

.  .  . 

10. 0 

6 

Laminar 

8.0 

5.0 

.  .  . 

12.1 

6.05 

Laminar 

8.0 

.  s 

3.1 

3.1 

Laminar 

8.0 

.5 

.  .  . 

5.6 

3.07 

Laminar 

8.0 

7.5 

.  .  . 

7.0 

2.8 

Laminar 

8.0 

15.0 

.  .  . 

6.0 

4.8 

Laminar 

8.0 

15.0 

.  .  . 

9.5 

Laminar 

8.0 

15.0 

... 

14.0 

Laminar 

37 

Tin 

16.0 

0 

m 

2.3 

1.7 

H 

■ 

Laminar 

X  o 

5 

3.5 

2.0 

Laminar 

lv  - 

7.5 

5.0 

2.95 

Laminar 

16.0 

10 

6.9 

4.7 

,  ■ 

ill 

Laminar 

16.0 

15 

Bfl 

27 

9.09 

HI 

Transition 

aNo  observed  separation. 
•^Detached  shock. 
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TRANSFER  COEFFICIENT 


PRESSURE  RATIO,  PM  /  Pr, 


PRESSURE  RATIO,  P?K/PFp 

FIG.  17.  Correlation  of  Peak  Heat-Transfer  Coefficient 
With  Peak  Pressure  for  Turbulent  Boundary  Layer. 
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4  6  3  10  20 

PRESSURE  RATIO,  P^/Pp* 


FIG.  18.  Correlation  of  Peak  Heat-Transfer  Coefficient 
With  Peak  Pressure  for  Transition. 
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FINAL  PRESSURE 

The  theoretical  inviscid  pressure-rise  across  a  shock  interaction, 
termed  the  final  pressure,  is  a  useful  correlating  parameter,  for  it 
describes  the  strength  of  the  interaction  and  requires  that  only  the 
upstream  Mach  number  and  the  configuration  be  known.  As  shown  in  Fig.  14, 
the  final  pressure  is  also  a  good  approximation  of  the  peak  pressure. 
Correlations  of  peak  heat-transfer  coefficient  with  final  pressure  for 
laminar,  turbulent,  and  transitional  boundary  layers  are  shown  in  Fig.  19, 
20,  and  21,  respectively.  Final  pressures  were  calculated  as  described 
in  the  Appendix.  Fin-shock  data  (Ref.  36  and  37)  are  not  included  in 
the  final  pressure  correlations  because  leading-edge  effects  on  the  flat 
plate  caused  the  Mach  number  upstream  of  the  interaction  to  differ  from 
the  Mach  number  generating  the  fin  shock. 

Because  peak  and  final  pressure  are  similar,  these  correlations 
show  the  same  trends  as  those  of  the  peak  pressure  correlations,  with 
the  turbulent  correlations  bring  considerably  better  than  the  laminar 
correlations. 
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FIG.  19.  Correlation  of  Peak  Heat-Transfer  Coefficient  With 
Theoretical  Final  Pressure  for  laminar  Boundary  Layer. 
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FIG.  20.  Correlation  of  Peak  Heat-Transfer  Coefficient  With 
Theoretical  Final  Pressure  for  Turbulent  Boundary  Layer. 
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FIG.  21.  Correlation  of  Peak  Heat-Transfer  Coefficient 
With  Theoretical  Final  Pressure  for  Transition. 


The 

pressure 


following  empirical  expressions  were  obtained  from  the  final 
correlations : 


Laminar: 


Turbulent : 


Transition: 


(based  on 
limited  data 


The  turbulent  expression  is  identical  to  that  obtained  in  the  peak 
pressure  correlations. 
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PLATEAU  PRESSURE 

It  is  suggested  in  Ref.  36  that  the  peak  heat-transfer  in  separated 
interactions  is  related  to  the  rise  from  plateau  pressure  to  final  pres¬ 
sure.  Correlations  of  peak  heat-transfer  coefficient  with  the  quantities 
Pf/^PL  an^  (?F  ~  PpL)/pFP  are  made  for  the  cases  in  which  separation  is 
known  to  exist.  Plateau  pressures  are  calculated  from  the  equations 
given  in  Section  4.  Laminar  data  correlated  best  with  the  quantity 
PF/PPL,  whereas  the  turbulent  data  were  best  described  by  (Pp  -  Pp^J/Ppp. 
Plots  of  these  two  correlations  appear  in  Fig.  22  and  23,  respectively. 
Plateau  pressure  correlations  were  not  made  with  the  transition  data. 
Correlating  peak  heat-transfer  with  plateau  pressure  provided  no  improve¬ 
ment  over  the  final  pressure  correlations.  As  in  the  previous  cases, 
the  laminar  data  exhibited  considerably  more  scatter  than  did  the  tur¬ 
bulent  data.  The  following  expression  relating  peak  heat-transfer  to 
plateau  pressure  in  turbulent  interactions  represents  a  least-squares 
fit  of  the  data  shown  in  Fig.  23.  An  expression  was  not  derived  for 
the  laminar  case  because  of  the  wide  scatter  in  the  data. 


REYNOLDS  NUMBER  AND  MACH  NUMBER 

In  ar.  attempt  to  obtain  better  peak  heat-transfer  correlations, 
additional  boundary-layer  parameters  were  introduced.  Parameters  used 
were  Reynolds  number  based  on  the  reference  distance  (Fig.  10),  Reynolds 
number  based  on  the  undisturbed  boundary-layer  thickness,  and  Mach 
number  upstream  of  the  interaction.  The  parameters  were  evaluated  as 
described  below. 

Reynolds  Number  Based  on  Distance 

Distance  Reynolds  numbers,  ReXR£,p,  were  based  on  free-stream  con¬ 
ditions  and  the  distance  from  the  leading  edge  of  the  model  to  the 
reference  point,  as  defined  in  Fig.  10.  As  laminar  and  turbulent 
Stanton  numbers  are  inversely  v'roportional  to  the  one-half  and  one-fifth 
powers,  respectively,  of  Reynolds  number,  these  values  were  used  in  the 
correlations. 
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FIG.  23.  Correlation  of  Peak  Heat-Transfer  Coefficient  With 
Final  and  Plateau  Pressures  for  Turbulent  Boundary  Layer. 


Reynolds  Number  Based  on  Boundary-Layer  Thickness 

Because  heat-transfer  coefficient  is  dependent  on  boundary-layer 
thickness,  it  was  felt  that  a  Reynolds  number  based  on  boundary-layer 
thickness  might  provide  a  better  correlation  than  would  distance 
Reynolds  number.  Reynolds  number  was  evaluated  by  using  the  hydro- 
dynamic  boundary-layer  thickness  at  the  reference  point.  The  boundary- 
layer  thickness  was  calculated  from  the  following  expressions  (Ref.  42): 


Laminar: 


Turbulent:  6 


°’376  XREF 
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The  boundary-layer-thickness  Reynolds  number,  Reg,  was  obtained  as 
follows : 


Re6 


6 


R  e6  = 


(- 


ZJW) 


Re6  = 


K 


(- 


V  x 


REF 


Re<5  = 


K 


( 


Re 

XREF 


Inserting  the  appropriate  values  for  the  constants,  a  and  K,  gives 


Laminar:  Re^ 


Turbulent:  Re 


As  in  the  case  of  the  distance  Reynolds  number  correlations,  the 
one-half  and  one-fifth  powers  of  boundary-layer-thickness  Reynolds 
number  were  used  in  the  laminar  and  turbulent  correlations,  respectively. 

To  attempt  peak  heat-transfer  correlations  with  the  change  in 
boundary- layer  thickness  across  the  interaction  region  would  have  been 
de.sirable.  However,  this  information  is  not  included  in  the  experimental 
data,  and  it  can  not  be  readily  obtained  analytically. 
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Upstream  Mach  Number 

In  most  of  the  experimental  data,  the  local  Mach  number  upstream 
of  the  interaction,  Mj ,  is  either  given  or  is  essentially  equal  to  the 
free-stream  Mach  number.  In  Ref.  30,  Popinski  calculated  the  local 
Mach  number  by  assuming  an  isentropic  compression  through  a  weak 
leading-edge  disturbance  from  the  free-stream  pressure  to  the  static 
pressure  on  the  plate  surface.  This  procedure  is  described  in  detail 
in  the  Appendix. 

The  correlations  attempted  by  using  Reynolds  number  and  Mach 
number  and  the  results  obtained  for  the  laminar  and  the  turbulent 
interactions  follow. 

Laminar  Interactions 


A  peak  laminar  heat-transfer  correlation,  of  the  form  suggested  by 
Popinski  (Ref.  30),  is  presented  with  the  use  of  the  parameters: 


hPK 


Re 


1/2  m3.7 
XREF  1 


PFP 


Popinski 's  correlating  parameters  differ  slightly  from  these  in  that 
he  used  Re,*,^  rather  than  Rex^Ep,  and  M^,  rather  than  M]_.  The  difference 

in  the  two  Reynolds  numbers  was  small,  and  it  was  felt  that  use  of  the 
Mach  number  upstream  of  the  interaction,  rather  than  the  free-stream 
Mach  number,  would  make  the  correlation  independent  of  leading-edge 
effects.  This  correlation,  shown  in  Fig.  24,  was  unsuccessful. 

Figure  25  shows  the  results  of  correlating  peak  heat-transfer  with 
final  pressure  and  Reynolds  number  based  on  the  undisturbed  boundary- 
layer  thickness.  This  correlation 


1/2 

5 


versus 


was  also  unsuccessful. 
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FIG.  25.  Correlation  of  Peak  Heat-Transfer  Coefficient  With 
Final  Pressure  and  Reynolds  Number  Based  on  Boundary  Layer 
Thickness  for  Laminar  Boundary  Layer. 


The  most  promising  of  laminar  correlations  vere  obtained  by  in¬ 
troducing  the  quantity  M^.  In  Fig.  26,  pea'/,  heac-transfer  coefficient 
was  correlated  with  and  final  pressure  i '  the  form 


2 

1 


versus 


In  Fig.  27,  the  Reynolds  number  based  or;  boundary-layer  thickness  is 
included  to  give  a  correlation  of  the  form 


Res1/2  M2 
hFP  S  1 


versus 


P 


FP 
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FIG.  26.  Correlation  of  Peak  Heat-Transfer  Coefficient  With 
Final  Pressure  and  Mach  Number  for  Laminar  Boundary  Layer. 


PRESSURE  RATIO,  Pr/?rr 

FIG.  27.  Correlation  of  Peak  Heat-Transfer  Coefficient  With 
Final  Pressure,  Mach  Number,  and  Reynolds  Number  Based  on 
Boundary-Layer  Thickness  for  Laminar  Boundary  Layer. 
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These  relations,  particularly  the  latter,  correlate  the  few  data 
points  other  than  those  of  Popinski  extremely  well.  The  Popinski  data 
appear  to  follow  the  same  slope  but  are  translated  upward,  indicating 
higher  heating  rates.  The  slope  of  these  data  indicates  a  relation  of 
the  form 


Turbulent  Interactions 

The  correlations  of  peak  turbulent  heat-transfer  coefficients  with 
peak  pressure  and  final  pressure  (Fig.  17  and  20,  respectively)  suggest 
that  peak  heat-transfer  increases  with  increasing  Mach  number.  The 
correlation  of  the  form 


versus 


was  attempted  to  account  for  this  effect.  The  results,  shown  in  Fig.  28, 
indicate  peak  heat-transfer  to  be  a  weaker  function  of  Mach  number  than 
that  attempted,  particularly  at  the  higher  Mach  numbers. 

Correlations  of  peak  turbulent  heat-transfer  with  Reynolds  number 
based  on  both  reference  distance  and  boundary-layer  thickness  provide 
a  slight  improvement  over  the  correlations  with  final  pressure  alone. 
Results  of  these  correlations  (Fig.  29  and  30)  were  used  to  develop  the 
following  two  relations : 
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FIG.  28.  Correlation  of  Peak  Heat-Transfer  Coefficient  With 
Final  Pressure  and  Mach  Number  for  Turbulent  Bourdary  Layer. 


PRESSURE  RATIO,  PF/PfF 

FIG.  29.  Correlation  of  Peak  Heat-Transfer  Coefficient  With 
Final  Pressure  and  Reynolds  Number  for  Turbulent  Boundary 
Layer. 
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Section  6 
RESULTS 


Results  of  the  peak  heat-transfer  correlations  indicate  the 
following  trends: 

1,  The  increase  in  heat-transfer  coefficient  caused  by  shock 
wave-boundary  layer  interactions  is  greater  in  laminar  interactions 
than  that  in  turbulent  interactions. 

2.  When  transition  occurs  in  the  interaction  region,  the  increase 
in  heat-transfer  coefficient  is  even  greater.  This  is  to  be  expected 
since  in  addition  to  the  increased-heating  rate  caused  by  the  shock 
interaction,  the  increasod-heating  rate  associated  with  turbulent  flow 
also  exists. 

Although  the  results  of  the  correlations  are  not  so  successful  as 
had  been  hoped,  particularly  in  the  case  of  laminar  interactions,  they 
can  be.  used  to  provide  a  rapid  preliminary  estimate  of  the  peak  heat- 
transfer  coefficient  in  shock  wave-boundary  layer  interaction  regions. 

For  this  purpose,  the  expressions  relating  peak  heat-transfer  to  final 
pressure  are  the  most  useful  of  all  developed  because  of  their  simplicity. 
The  more  involved  correlations  provided  little  improvement  and,  in  some 
cases,  were  not  so  successful  as  the  correlations  with  fir.al  pressure 
alone . 

Possible  explanations  for  the  lack  of  success  achieved  in  the  peak 
heat-transfer  correlations  follow. 


EXPERIMENTAL  TECHNIQUES 

Although  all  of  the  data  used  in  the  correlations  came  from  wind 
tunnel  investigations,  in  which  test  conditions  may  be  carefully 
monitored,  several  variations  in  experimental  techniques  were  involved. 
Some  of  the  studies  utilized  separate  models  for  measuring  pressure  and 
heat-transfer  distributions  so  that  these  parameters  could  not  be  mea¬ 
sured  simultaneously.  Correlations  between  pressure  and  hea  transfer 
were  thus  affected  by  possible  variations  in  flow  condition 

Although  most  of  the  studies  did  not  state  the  overall  ac.  racy  of 
the  data,  several  did  mention  difficulties  encountered  during  testing. 
Experimental  accuracies  of  15  and  20%  were  given,  and  it  is  doubtful  if 
any  of  the  data  are  much  more  accurate  than  this.  This  could  account 
for  much  of  the  scatter  in  the  correlations. 
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Data  were  obtained  for  interactions  on  both  hot  and  cold  walls  and 
under  transient  and  steady-state  conditions.  The  difference  in  heat- 
transfer  measurement  techniques  and  the  various  assumptions  made  in 
reducing  the  data  to  local  heating  rates,  Stanton  numbers,  and  heat- 
transfer  coefficients  added  to  the  difficulty  of  comparing  the  results 
of  various  tests. 


MEASUREMENT  OF  PEAK  HEATING  RATE 

Because  large  gradients  in  heating  rate  exist  in  shock  wave-boundary 
layer  interaction  regions,  it  is  difficult  to  measure  peaks  in  heating- 
rate  distributions.  When  the  spacing  of  sensors  is  not  dense  enough,  it 
is  possible  to  miss  a  measurement  at  or  near  the  peak  value. 


METHOD  CF  PRESENTING  DATA 

Ideally,  all  of  the  heat-transfer  data  should  be  nondimens ionalized 
by  a  common  reference  value,  such  as  a  theoretical  heat-transfer  coef¬ 
ficient  calculated  by  a  particular  method.  Unfortunately,  because  the 
published  data  were  presented  ii  a  variety  of  ways,  this  was  not  possible. 
It  was  desired  to  reference  all  peak  heat-transfer  coefficients  to  the 
value  that  would  exist  at  the  reference  point  (Fig.  10)  if  the  flow  were 
undisturbed.  Reference  values  included  theoretical  and  experimental 
heat-transfer  coefficients  for  undisturbed  flow  over  flat  plates,  and 
experimental  values  upstream  of  the  interaction.  In  most  of  the  cases 
in  which  the  letter  value  was  used,  the  difference  between  the  reference 
point  and  the  beginning  of  the  interaction  region  was  small,  compared  tc 
the  distance  from  the  beginning  of  the  boundary  layer.  Consequently, 
the  difference  in  heat-transfer  coefficient  at  the  two  locations  was 
also  small. 

In  some  cases,  it  was  difficult  to  accurately  determine  the  peak 
values  of  the  pressure  and  heat-transfer  distributions  from  the  graphical 
forms  used  in  presenting  the  data. 


STATE  OF  THE  BOUNDARY  LAYER 

A  major  cause  for  the  poor  correlations  obtained  with  the  laminar 
data  might  be  the  possibility  of  undetected  transition  in  the  inter¬ 
action  region.  This  could  explain  why  the  data  of  Popinski,  in  Fig.  26 
and  27,  are  translated  upward  from  the  rest  of  the  data  points,  which 
correlate  quite  well. 

Another  cause  of  scatter  in  the  correlations  might  be  the  nature 
of  separated  flow,  which  makes  if.  difficult  to  repeat  results  on 
successive  test  runs  let  alone  to  compare  results  obtained  at  different 
facilities. 
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Section  7 

CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FURTHER  STUDY 


CONCLUSIONS 

1.  Heating  rates  in  shock  wave-boundary  layer  interaction  regions 
are  significantly  higher  than  ♦'he  normally  predicted  values.  These 
increase j  in  heating  rate  are  highest  for  cases  in  which  transition 
occurs  in  the  interaction  region.  Greater  heating  rate  increases  occur 
in  laminar  interac cions  than  in  turbulent  interactions. 

2.  •  Jficient  experimental  data  are  not  available  to  allow  the 
development  of  rapid  empirical  means  for  predicting  the  heat-transfer 
distribution  in  interaction  regions.  Correlations  are  made  between 
peak  experimental  heating  rates  in  interaction  regions  and  various  flow- 
field  parameters  with  limited  success. 

3.  Peak  heat-transfer  in  turbulent  interaction-regions  may  be 
estimated  by  the  following  relation  that  involves  the  inviscid  pressure 
rise  across  the  interaction: 


Peak  heat-transfer  data  for  laminar  interactions  do  not  correlate 
nearly  so  well  as  the  turbulent  data;  the  following  expression  is  sug¬ 
gested  for  use  in  providing  qualitative  estimates  only: 


Sufficient  experimental  data  arc  not  available  for  transitional 
interactions  tc  al’^w  development  of  empirical  peak  heat-transfer 
prediction  techniques. 


RECOW^  DAT  IONS  FOR  FURTHER  STUDY 

1.  The  geometry  of  shock  wave-boundary  layer  interaction  regions 
(interaction  and  separation  lengths)  should  be  further  investigated  to 
allow  the  deve.  pment  of  general  expressions  for -predicting  heat-transfer 
distributions. 
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2.  More  experimental  data  are  required-part-1  cularly  for  laminar 
and  transitional  interactions— so  that  more  valid  expressions  for 
predicting  peak  heating  rate  may  be  developed.  Extreme  care  should 
be  taken  in  conducting  laminar  investigations  to  determine  whether  or 
pot  transition  occurs  in  the  interaction  region. 

3.  Investigation  of  the  boundary-layer  thickness  through  the 
interaction  region  would  provide  a  useful  parameter  for  peak  heat- 
transfer  correlations. 

4.  Studies  should  be  made  to  investigate  heat- transfer  distributions 
immediately  downstream  of  interaction  regions.  This  would  be  particularly 
useful  for  shock-impingement  interactions. 
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Appendix 

CALCULATION  OF  THEORETICAL  INVISCID  PRESSURE 
RATIO  ACROSS  SHOCK  INTERACT70N 

The  theoretical  inviscid  pressure  ratio  across  a  shock  interaction 
nay  be  calculated  by  using  the  applicable  oblique  (Ref.  46)  or  conical 
(Ref.  47)  shock  relations.  These  relations  require  that  only  the  Mach 
number  upstream  of  the  interaction,  M^,  and  the  flow  deflection  angle, 
9,  be  known.  In  the  case  of  a  compression-corner  interaction  with 
separation,  an  improved  estimate  of  the  pressure  ratio  may  be  made  by 
including  Reynolds  number  effects. 

If  the  upstream  Mach  number  is  not  known,  it  may  be  calculated 
ftom  the  upstream  pressure,  free-stream  pressure,  and  free-stream  Mach 
number.  The  following  procedure  assumes  an  isentropic  compression  of  a 
perfect  gas  (y  =  1.4)  through  a  weak  leading-edge  disturbance: 


PT  =  PT  (*sentroPic  compression  from  P^  to  P^) 
1  oo 


/  2  \"7/2 

00  \i  +  H2J5J 
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The  methods  used  to  calculate  inviscid  pressure  ratios,  expressed 
as  Pp/Ppp  in  the  correlations  given  in  Section  5,  are  described  for  the 
various  configurations  in  the  following  text.  Nomenclature  refers  to 
Fig.  10. 


SHOCK  IMPINGEMENT  ON  FLAT  PLATE 

The  pressure  ratio,  Pf/**FP»  is  calculated  using  the  oblique  shock 
rel<»  '  -sns: 


6X  -  f(Mx  ,  0) 


and  is  obtained  by  solving  the  following  cubic  equation: 


where 


sin^  +  b  sin^  8^  +  c  sin2  8^  +  d  *  0 


(1) 


2 

M^  +  2 


-  y  sin  9 


2  Mj  +  1 


(Y  +  1)  Y  -  1 
4  M? 


sin2  6 


d  -  - 


2 

cos  8 
M4 


P2  7  M2  sin2  02  -  1 


=  s. 


(2) 
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Cl+6  1/2  3 

2  7  sin2  (B1  -  9) 

S2  »  f (M2  ,  6)  (Ref.  Eq.  1) 

P,.  7  M2  sin2  B0  -  1 


SHOCK  IMPINGEMENT  ON  CYLINDER 

Calculations  for  conditions  in  Region  2  (8]_,  P2/P1,  M2)  ar*»  the 
same  as  those  for  impingement  on  a  flat  plate. 

Conditions  in  Region  3  (63,  P3/P2>  M3)  are  obtained  from  empirical 
conical  shock  relations  (Ref.  47). 
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If  0  -  sin"*  (G2/G3)^  in  the  above  equations,  the  shock  is  detached, 
and  a  solution  is  not  readily  obtainable. 


*[' 


F2  +  ¥  sin  0 


-  |[f2  -  si„2  e] 


where 
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1 


Y  +  7 
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3 


,  1  /Y  +  7 

2  \Y  +  1 


f3s 

p2 


+  1 


(8) 


(9) 


(10) 


The  Mach  number  on  the  surface  of  the  cylinder,  M3C,  is  obtained 

by  assuming  an  isentropic  compression  from  the  pressure  and  the  Mach 
number  immediately  behind  the  reflected  conical  shock,  P3  and  M3,  to 
the  pressure  on  the  surface  of  the  cylinder,  P3  .  From  oblique  shock 
relations  P3  and  M3  are  calculated: 


k 

i 

1 
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(11) 


WEDGE  COMPRESSION  CORNER 
Attached  Boundary  Layer 

The  calculation  of  8].,  p2/,pl»  ant*  M2  same  as  that  for  the 

case  of  shock  impingement  on  a  flat  plate: 


Separated  Boundary  Layer 

In  the  case  of  separated  compression-corner  flow,  the  overall 
pressure  rise  is  assumed  to  take  place  by  means  of  a  double  compression. 
The  pressure  behind  the  separation  shock  wave,  P2»  is  assumed  to  be 
identical  with  the  plateau  pressure,  pPL»  and  is  obtained  from  the 
expressions 


CpPHAM  *  l-‘  KrEF )'1M  (M1  -  4  (12) 

CPPLTHRB  ’  1'7  (ReXREF)'1/10  (MX  '  4  ^  ‘13> 

"  (CPPL  H1  l)  +  1  ’  h  <“> 


63 


NWC  TP  4485 


Equations  12  and  13  are  taken  from  Ref.  43.  The  separation  angle,  9', 
is  calculated  as  being  the  deflection  angle  required  to  create  the  rise 
to  the  plateau  pressure. 


6  £  +  1 


,1/2 


..  _  -1  I  "1 


0'  =  sin 


7  Mt 


M. 


Mi  (6  h  +  l)  -  5  f5i  -  0 


ll/2 


(?i  +  6) 


sin 


k  +  6  \1/2 

K  - e,)  ■  M 


e"  =  e  -  e’ 


(15) 


(16) 


(17) 

(18) 


B2  =  f (M2  ,  e")  (Ref .  Eq.  1) 

P3  7  M2  Sin2  e2  -  1 


^2  +  6 
7  C2  sin2  (82  -  9") 


2 

1/2 
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CONICAL  FLARE  uOMPRESSION-CORNER 


Attached  Boundary  Layer 


Conditions  in  Region  2  (82,  p2/pl»  p2s*  M2>  anc*  M2g^  are  calculated 

by  using  the  same  equations  as  those  used  for  shock  reflection  from  a 
cylinder. 


6,  =  f (Mj  ,  0) 


(Ref.  Eq.  7) 


=  f(K’  ,  0) 

1 


(Ref.  Eq.  8  and  9) 


*r  =  f(Kl  ’  Bi) 
1 


(Ref.  Eq.  2) 


M2  =  f(Bx  ,  0) 


(Ref.  Eq.  3) 


(Ref.  Eq.  11) 


Separated  Boundary  Layer 


The  Ideal  method  of  calculating  the  separation  angle,  0’,  would  be 
to  assume  that  p2g/pi  is  equal  to  Ppl/Pj.,  to  calculate  ppl/p1  from 

Eq .  12  or  13,  and  to  calculate  0'  from  Eq.  8.  Because  there  is  difficulty 
in  solving  Eq.  8  for  0',  the  overall  pressure  rise  may  be  approximated  by 
assuming  that  0'  =  8"  =  0/2.  The  procedure  for  calculating  Pp/Ppp  then 


becomes 
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t 


Cp3s  ' 


f(Mi  , 

6') 

(Ref.  Eq.  7) 

f(MX  , 

0') 

(Ref.  Eq.  8) 

f(Cp2S 

.  Ml) 

(Ref.  Eq.  9) 

f(Ml  , 

Si) 

(Ref.  Eq,  2) 

Hh  , 

6') 

(Ref.  Eq.  3) 

/P2s 
\P2  ’ 

M2j 

(Ref.  Eq.  11) 

f(Mps  ' 

0") 

(Ref.  Eq.  7) 

f(H2s  , 

0"j 

(Ref.  Eq.  8) 

f(Cp3S 
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(Ref.  Eq.  11) 


(20) 


FIN  SHOCK 

The  calculation  of  02>  ^2^1*  an<^  ^2  t^ie  same  as  that  for  shock 
impingement  on  a  flat  plate. 


P 


F 


P 


FP 


METHODS  OF  SOLUTION 

Graphical  solutions  are  available  for  all  of  the  previously  mentioned 
parameters  (Ref.  46).  The  large  number  of  data  points  involved  in  the 
correlations,  however,  made  graphic-',  solutions  impractical.  For  this 
reason,  a  FORTRAN  IV  computer  progr  -  was  written  for  the  UNIVAC  1108 
computer  to  solve  the  inviscid  flow  equations. 
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